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CHAPTER  1 
INTRODUCTION 


Water  management  is  a  crucial  element  in  directing 
future  land  development  across  the  united  States.  Measures 
are  often  taken  to  minimize  flood  hazards  while  at  the  same 
time  maintaining  an  adequate  water  supply  for  consumption, 
irrigation,  sewage  treatment  and  aquifer  regeneration 
(Stockton,  1975) .  water  resource  management  activities 
directed  at  regulating  stream  runoff  must  address  a  wide 
variety  of  considerations,  including  the  frequency  and 
magnitude  of  floods.  The  confidence  of  flood  frequency 
analysis  is  dependent  on  the  existence  and  length  of  gaged 
records  (Kochel  and  Baker,  1982;  Martens,  1992) .  A  lack  of 
long- duration  hydrologic  records  for  United  States  streams 
has  limited  the  confidence  of  flood  frequency  estimates,  and 
created  am  opportunity  for  data  contributions  by  proxy 
methods . 

The  science  of  tree- ring  amalysis,  or  dendrochronology, 
cam  provide  a  valuable  source  of  paleoflood  information 
(Martens,  1992;  McCord,  1990;  Sigafoos,  1964).  The  year  and, 
often,  magnitude  of  significant  floods  can  be  determined  from 
the  amnual  ring  by  considering  the  geomorphic  situation  and 
biological  response  of  trees  to  flooding.  The  growth  of  tree 
rings  primarily  depends  on  soil  nutrients  (nitrogen, 
phosphorous,  potassium,  calcium,  magnesium,  and  sulfur),  and 
climatic  factors  (mainly  temperature  and  soil 
moisture) (Kramer  and  Kozlowski,  1960).  By  calibrating  ring 
growth  to  recorded  weather  and  discharge  data,  tree -ring 
information  can  be  used  to  reconstruct  yearly  variations  in 
climate.  These  relationships  can  be  applied  through  the  life 
of  the  tree,  thus  extending  recent  climatic  and  hydrologic 
records  back  in  time  (Fritts,  1976). 
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l.A.  Purpose  and  scope 


The  purpose  of  this  study  is  to  demonstrate  how  tree¬ 
ring  analysis  along  rivers  can  provide  valuable  pre -historic 
flood  frequency  and  magnitude  information.  The  Big  Muddy 
River  of  southern  Illinois  w*  3  selected  as  a  candidate  for 
paleoflood  research.  Justification  for  its  selection,  and 
specific  objectives  and  hypotheses  follow. 

l^A.l _ Justification  for  research 

The  Big  Muddy  River,  a  tributary  of  the  Mississippi 
River,  lies  in  a  region  which  is  historically  prone  to 
flooding.  Major  floods  since  the  turn  of  the  century  have 
resulted  in  significant  loss  of  property  and  life,  and 
prompted  numerous  flood  control  projects.  The  scope  of  these 
projects  were  based  on  flood  information  from  gaged  records, 
which  did  not  generally  begin  in  the  area  in  earnest  until 
after  1910.  To  base  flood  magnitude  and  frequency  estimates 
on  such  a  limited  observed  record  is  insufficient,  since  the 
variance  of  the  gaged  period  may  not  be  representative  of  the 
variance  over  longer  time  periods.  Long  period  variance  is 
often  the  best  predictor  for  future  events,  since  it  has 
increased  potential  to  encompass  extremes  and  minor 
fluctuations  which  may  someday  be  repeated. 

To  date,  there  has  been  little  effort  to  reconstruct 
long  period  variance  of  discharge  on  the  Big  Muddy  River,  or 
in  the  upper  Mississippi  River  valley.  Such  research  has 
potential  economic  value,  since  increasing  land  development 
is  forcing  drainage  districts  to  assess  their  ability  to 
handle  high  flow  events,  while  maintaining  adequate  water 
levels  when  precipitation  is  low.  Large  tracts  of 
residential,  commercial,  and  agricultural  land  lie  on 
floodplains  and  are  vulnerable  to  floods,  so  it  is  important 
to  improve  our  understanding  of  flood  potential.  But  since 
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resources  are  limited,  the  costs  and  benefits  of  flood 
control  projects  must  be  carefully  weighed.  The  most 
important  elements  in  project  design  are  flood  magnitude  and 
frequency,  which  together  can  be  used  to  calculate  risks  that 
future  flood  events  will  exceed  protection  limits.  An 
extension  of  the  hydrologic  record  using  tree  rings  can 
provide  more  complete  information  on  flood  risks,  and 
therefore  better  direction  for  measured  flood  control. 

l.A.2  Objective 

The  objective  of  this  research  is  to  reconstruct  past 
river  discharge  and  flood  events  on  the  Big  Muddy  River  using 
dendrochronology.  Suitable  trees  were  located  in  the  Shawnee 
National  Forest,  and  climatic  and  hydrologic  data  for  the 
past  60  years  were  available  for  initial  calibration. 

The  following  questions  are  addressed:  How  does  tree 
growth  at  various  levels  on  and  above  the  floodplain  reflect 
flood  events?  What  do  these  relationships  infer  about  past 
river  behavior?  And  finally,  what  applications  are  there  to 
other  developing  areas  in  Illinois  or  the  Upper  Mississippi 
River  valley? 

JL»A.3 _ Hypotheses 

The  general  working  hypothesis  for  this  research  is  that 
flow  history  can  be  reconstructed  from  the  tree-ring  record 
by  developing  a  model  which  correlates  recent  tree- ring 
growth  with  recent  hydrologic  data.  The  sampling  and 
measurement  of  tree  cores,  statistical  manipulation  of  ring 
widths,  and  model  development  will  be  explained  later. 

Hypotheses  specific  to  the  study  are  that  tree- ring 
analysis  on  floodplains  and  terraces  enables: 
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a.  Determination  of  a  prehistoric  flood  chronology:  Years 
in  which  flooding  occurred  are  identified  for  the  period  of 
the  tree- ring  chronology. 

b.  Estimation  of  flood  heights:  Heights  of  flood  scars  on 
trees  above  the  floodplain,  elevation  of  flood- ringed  trees 
above  the  floodplain,  and  comparison  of  tree- ring  indices 
from  different  levels  on  and  above  a  floodplain  serve  as 
rough  estimates  of  flood  heights  for  the  field  site  river, 
assuming  that  significant  floodplain  degradation  or 
aggradation  did  not  occur  during  previous  decades. 

c.  Reconstruction  of  mean  monthly  discharge:  Correlation 
and  regression  of  tree- ring  growth  on  gaged  discharge  data 
provides  a  mathematical  means  of  estimating  mean  monthly 
discharge  from  ring  width. 

1.A.4 _ Organization  of  the  thesis 

Chapter  1  presents  the  study  hypotheses,  and  provides 
background  information  on  how  paleoflood  information  is 
stored  in  tree  rings.  Current  methods  and  applications  of 
dendrochronology  are  discussed. 

Chapter  2  introduces  the  study  site  where  flood  research 
was  conducted.  A  geomorphic  analysis  and  recent  flood 
history  of  the  area  is  included. 

Chapter  3  describes  methods  of  data  acquisition  and 
analysis.  Tree  core  collection  is  presented  along  with 
sources  and  quality  of  hydrologic  and  climatic  data.  Model 
development  and  related  statistics  are  outlined. 

Chapters  4  and  5  discuss  the  results  and  conclusions  of 
the  flood  reconstruction  project,  as  well  as  the  ioqpact  of 
such  information  on  future  planning. 

Before  proceeding  with  any  scientific  research,  it  is 
important  to  have  an  understanding  of  the  basic  concepts  and 
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methods  in  the  area  of  study.  Therefore,  the  remainder  of 
this  introductory  chapter  is  devoted  to  outlining  important 
points  in  dendrochronology,  modern  flood  analysis,  and 
paleoflood  analysis. 

l.B  Concepts  in  dendrochronology 

A  review  of  prevalent  themes  in  dendrochronology 
outlines  the  current  framework  for  climatic  reconstruction. 
Concepts  in  dendrochronology  related  to  paleohydrology  are 
discussed  to  provide  direction  for  future  dendrohydrologic 
research  on  Illinois  river  systems. 

l.B.l  Tree  response  to  climatic  inputs 

Tree- ring  growth  represents  a  year-long  integration  of 
its  environment  (Fritts,  1965) .  Although  total  seasonal 
growth  is  the  result  of  many  interacting  factors,  only  one 
environmental  factor  dominates  in  limiting  growth,  i.e. , 
available  soil  moisture  during  the  growing  season  (Stokes  and 
Smiley,  1968;  LaMarche,  1974;  Phipps,  1972;  Stockton,  1975). 
Wider  rings  generally  correspond  to  a  moist,  cool  year;  while 
narrow  rings  depict  dry  and/or  warm  conditions  (Fritts, 

1965) . 

Ring  widths  are  determined  by  a  yearly  flush  of  growth 
that  begins  in  spring  and  ends  in  summer  (Figure  1) .  During 
the  early  part  of  the  growing  season,  large  "earlywood”  cells 
are  produced.  Toward  the  middle  and  end  of  the  season, 
growth  slows,  and  smaller,  thin-walled  NlatewoodN  cells  are 
formed.  The  boundary  between  annual  rings  is  marked  by  the 
abrupt  shift  from  latewood  to  earlywood  cells.  The  width  and 
character  of  the  annual  rings  is  dependent  not  only  on 
climatic  inputs,  but  site  and  species  characteristics  as 
well. 
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Earlywood  Latewood 


Annual  growth  ring 


Annual  ring  growth  (from  Stokes  and  Smiley,  1968) 
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Site  selection  in  paleoflood  research  aims  to  sample 
trees  whose  growth  is  sensitive  to  changes  in  hydrologic 
parameters.  Trees  on  well  drained  sites  (such  as  slopes  or 
terraces)  are  more  sensitive  to  changes  in  precipitation 
since  they  are  shallow- rooted  and  readily  affected  by  a  rapid 
change  in  soil  moisture  supply  (Stockton,  1975;  Smith  and 
Stockton,  1981) .  These  trees  tend  to  exhibit  a  varying,  or 
"sensitive"  ring- width  series  (Figure  2) .  Poorly  drained 
sites  (such  as  floodplains)  usually  realize  optimum 
conditions  for  growth,  and  a  "complacent"  ring-width  series, 
which  varies  little  from  year  to  year.  Complacency  in 
floodplain  trees  can  be  interrupted  if  flooding  saturates  the 
soil  and  limits  the  delivery  of  oxygen  to  the  roots  (Fritts, 
1965b) .  Lack  of  oxygen  in  the  soil  leads  to  accumulation  of 
carbon- dioxide  and  toxic  metabolic  products  that  can  inhibit 
growth  of  the  tree  or  cause  injury  to  living  cells  (Gill, 
1970).  The  degree  of  flood- related  impact  on  a  tree's  growth 
is  highly  dependent  upon  tree  species  and  timing  and  duration 
of  flooding  (Gill,  1970;  Hosner  and  Boyce,  1962). 

1.BL.3 _ The  importance.  ol .  species  selection 

Tree  species  selection  is  important  since  various 
species  do  not  respond  to  flooding  in  the  same  way.  Also, 
bottomland  trees  have  different  tolerances  to  saturated  soils 
(Hosner  and  Boyce,  1962;  Phipps,  1972).  Seme  species  thrive 
in  saturated  conditions  while  others  suffer  limited  growth  or 
even  death.  In  addition,  tree  selection  for  coring  and  study 
must  be  grouped  by  species  to  eliminate  species -related 
growth  variations  (Fritts,  1976) . 

For  most  deciduous  species,  flooding  affects  ring  growth 
only  when  it  occurs  during  the  growing  season  (Gill,  1970; 
McAlpine,  1961) .  This  is  mainly  because  oxygen  demand  by  the 
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Figure  2.  Site  influence  on  tree  sensitivity  (from  Stokes 
and  Smiley,  1968) . 


roots  is  minimal  in  winter  (Gill;  1970) .  Dormant  season 
flooding,  however,  may  lead  to  increased  soil  moisture  for 
the  start  of  the  following  growing  season.  Duration  of 
flooding  during  the  growing  season  is  very  important,  as 
tolerances  vary  dramatically  with  species  and  tree  age. 
Seedlings  and  young  trees  are  most  sensitive  to  inundation, 
since  their  root  systems  are  less  developed  (Gill,  1970; 
Hosner  and  Boyce,  1962) . 

In  addition  to  understanding  the  biological  response  of 
trees  to  flooding,  it  is  important  to  conprehend  magnitude 
and  frequency  estimation. 

l.C  Magnitude  and  frequency  of  floods 

A  key  c exponent  to  floodplain  management  is  flood 
frequency  estimation,  which  is  the  probability  of  occurrence 
for  a  flood  to  exceed  a  given  magnitude  (Knighton,  1984) . 
Recurrence  interval,  the  reciprocal  of  frequency,  is  the  mean 
number  of  years  separating  flood  events. 

The  conventional  method  of  flood  frequency  estimation 
uses  an  annual  maximum  flood  series,  in  which  the  highest 
peak  discharge  of  each  water  year  is  listed  in  an  array  that 
is  ranked  to  size  from  largest  to  smallest.  The  return 
period  for  a  particular  year's  peak  discharge  is  computed  as 

T=(n+l) /m 

where  T  is  return  period  (years) ,  n  is  number  of  years  of 
record,  and  m  is  rank  within  the  series  (Knighton,  1984) . 

The  return  period  estimate  is  limited  by  the  length  of 
the  flood  series  record,  with  a  maximum  return  period  of  n+l 
years  for  the  largest  (m=l)  flood  in  the  series.  Return 
periods  greater  than  n+l  years  can  be  estimated  by 
graphically  extending  a  logarithmic  flood  frequency  curve  to 
the  desired  probability  (Kite,  1977).  Uncertainty  is  added 
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by  assuming  that  the  form  of  the  extending  curve  remains  the 
same  for  rare  events.  In  addition,  these  return  period 
estimates  could  be  inaccurate  because  the  period  of  gaged 
discharge  may  not  be  representative  of  longer  period 
discharge.  The  gaged  period  may  be  shorter  than  return 
periods  of  particularly  destructive  floods,  or  coincide  with 
anomalous  climatic  conditions  (Kochel  and  Baker,  1982; 
Hirschboeck,  1987).  To  illustrate,  a  calculated  50  year 
flood  event  from  a  49  year-long  gaged  record  could  in  fact  be 
a  100  or  200  year  event  that  just  happened  to  occur  in  the 
gaged  period.  There  is  no  guarantee  that  the  expected  number 
of  "return"  years  will  separate  flood  events.  Ideally,  a 
time  series  should  be  long  enough  to  encompass  the  full  range 
of  discharge  fluctuations.  Realistically,  gaged  records  must 
be  extended  with  proxy  data  to  make  return  period  estimates 
with  greater  confidence. 

A  proven  source  of  proxy  flood  data  is  dendrochronology. 
Current  techniques  will  be  discussed  in  the  remainder  of  this 
chapter. 

l.D  Paleoflood  analysis  using  tree  rings 

Dendrochronology  can  potentially  provide  two  types  of 
hydrologic  information:  l)  years  of  extreme  hydrologic  events 
such  as  floods  or  droughts,  and  2)  a  paleoclimatic  record 
useful  for  extending  hydrologic  records  and  estimating  trends 
and  variability  of  precipitation  and  discharge  (Stockton, 

1975) .  A  number  of  analysis  techniques  exist  and  have  been 
proven  effective.  The  most  common  approaches  are  described 
below: 

1J3.1 _ Flood  rings 

Years  of  prolonged  flooding  in  riparian  areas  cam  be 
determined  with  the  discovery  of  "flood  rings".  Flood  rings 
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are  identified  by  abnormalities  in  the  wood  tissue  of  the 
flooded  season’s  growth  ring  (Gill,  1970;  Yanosky,  1983; 
Martens,  1992) .  The  abnormalities  appear  as  enlarged  cells, 
similar  to  earlywood  cells,  but  formed  in  a  non - ear lywood 
portion  of  the  growth  ring  (Figure  3) .  The  enlarged  cells 
form  when  l)  trees  produce  a  second  crop  of  leaves  during  the 
growing  season,  or  when  2)  rapid  shoot  elongation  occurs  and 
new  leaves  are  formed  at  stem  tips  (Yanosky,  1983;  Phipps, 
pers.  comm.,  1993). 

The  first  flood  ring  type  occurs  with  some  degree  of 
crown  inundation,  where  trees  are  partially  or  totally 
defoliated.  If  injury  occurs  early  enough  in  the  growing 
season,  growth  will  continue  with  formation  of  new  leaves. 
Hormonal  release  from  refoliation  often  causes  development  of 
large  vessels  in  the  growth  ring  (Yanosky,  1983) . 

The  second  type  of  flood  ring  is  also  related  to 
flooding,  though  non- injurious  to  a  tree.  Root  inundation 
for  a  few  days  may  result  in  rapid  growth  fueled  from  added 
moisture  (Phipps,  pers.  comm.,  1993).  Hormones  released  from 
new  leaf  growth  at  stem  tips  may  cause  a  band  of  large  cells 
(called  "white"  rings)  to  be  formed  in  the  growth  ring,  often 
as  an  extension  of  the  earlywood  zone  (Yanosky,  1983) . 

Flood  rings  can  be  discovered,  should  they  exist, 
through  extensive  tree  coring  of  floodplain  forests.  Cores 
from  about  10  or  more  flood-prone  trees  must  be  carefully 
examined  under  a  microscope  for  growth  abnormalities. 
Annotation  of  abnormal  growth  years  shows  flood  frequency, 
highlighting  common  years  of  flood  evidence  from  floodplain 
trees.  The  time  of  year  when  flooding  occurred  can  be 
determined  by  noting  the  intra-ring  position  of  the  flood 
ring  relative  to  the  total  ring  width,  which  gives  the 
approximate  time  of  flooding  during  the  growing  season 
(Yanosky,  1983) .  Comparison  of  tree  base  heights  on  uneven 
floodplains  may  give  information  on  flood  magnitude,  once 
elevations  of  affected  trees  are  identified  for  flood  years. 
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1.D.2  Flood  scars 

Years  of  flooding  can  also  be  determined  through  flood 
scar  analysis.  High  discharge  events  have  the  ability  to 
transport  large  quantities  of  debris  at  high  velocities. 

When  over-bank  flooding  occurs,  floodplain  trees  are 
vulnerable  to  physical  damage  on  the  upstream  side  by 
floating  logs  or  other  materials  rolling  on  the  bed.  Damage 
is  most  likely  to  occur  at  peak  discharge,  when  logs  are 
being  entrained  into  the  flow.  Trees  adjacent  to  the  channel 
are  at  highest  risk  since  logs  are  first  intercepted  there, 
and  velocity  slows  further  inland  (Gottesfeld  and  Johnson - 
Gottesfeld,  1990) .  The  most  common  types  of  tree  damage  are 
uprooting,  breaking,  shearing,  and  scarring  (Sigafoos,  1964) . 

A  scarring  blow  to  the  tree  stem  results  in  injury  to 
the  bark  and  vascular  cambium  (Figure  4) .  The  bark  and  dead 
tissue  are  naturally  sloughed  off  and  replaced  by  callus 
tissue,  which  begins  to  cover  the  exposed  wood  from  the 
margins  (McCord,  1990;  Sigafoos,  1964) .  Eventually  the 
scarred  area  is  healed  over,  and  couplet e  annual  rings  once 
again  encircle  the  stem. 

Tree  scars  can  be  dated  by  either  of  two  methods.  The 
first  approach  is  to  core  the  damaged  tree  and  count  the 
number  of  rings  in  from  the  bark  to  the  last  sound  ring 
before  scarring.  Dating  is  also  possible  after  coring 
through  the  scar,  marking  the  damaged  ring,  and  cross -dating 
that  ring  with  undamaged  sections  of  the  tree. 

Flood  scar  analysis  was  applied  by  Sigafoos  (1964)  to 
date  flood  events  on  the  Potomac  River  in  Virginia,  where 
evidence  from  28  floods  were  found  in  stems  of  ash  trees. 
Techniques  can  also  be  applied  to  remote,  poorly- gaged 
rivers.  Major  flood  events  and  channel  changes  over  the  past 
105  years  were  documented  for  the  Morice  River,  a  braided 
gravel-bed  river  in  British  Columbia,  using  flood  scars  and 
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Figure  4.  Illustrations  of  (a)  undamaged  tree  (b)  tree 
injured  by  physical  inpact  from  flooding,  (c,d)  formation  of 
callus  tissue,  (e)  healed  flood  scar  (after  Sigafoos,  1964) . 


forest  stand  establishment  dates  (Gottesfeld  and  Johnson- 
Gottesfeld,  1990) . 

Measurements  of  scar-heights  are  a  source  of  paleostage 
data,  which  enables  a  crude  estimate  of  past  discharge. 
Surveyed  channel  dimensions  and  scar-heights  in  Rattlesnake 
Canyon,  Arizona,  have  been  used  to  to  compute  paleo discharge 
for  a  number  of  flow  depths  (McCord,  1990) . 

1.D.3  Floodplain- upland  comparisons 

The  effects  of  discharge  on  tree- ring  growth  can  be 
isolated  by  analyzing  ratios  of  annual  floodplain  tree- ring 
growth  over  annual  upland  growth  (Wendland  and  watson- 
Stegner,  1983) .  This  technique  is  based  on  the  premise  that 
while  trees  from  both  sites  experience  the  same  climate, 
trees  located  on  flooded  soil  potentially  have  comparatively 
wider  or  narrower  rings  than  upland  trees  .  The  added  soil 
moisture  from  flooding  may  lead  to  a  flush  of  growth  and  an 
abnormally  wide  ring.  With  extended  duration  floods, 
prolonged  soil  saturation  during  or  just  prior  to  the  growing 
season  may  cause  a  lack  of  oxygen  transfer  to  the  roots  and 
inhibited  ring  growth  for  that  year  (Mitsch  and  Rust,  1984; 
Gill,  1970;  Hosner  and  Boyce,  1962) .  A  floodplain/upland 
ratio  near  1.0  indicates  no  flooding  for  that  growing  season. 
A  ratio  which  is  signif icantly  less  or  greater  than  1.0  shows 
that  one  or  more  growth-altering  floods  were  present  on  the 
floodplain  during  the  spring  or  simmer  of  that  year. 

This  type  of  study  requires  core  sampling  from  the 
floodplain  and  a  number  of  topographic  elevations  above  it. 
Tallying  years  with  ratios  less  than  one  enables  flood 
frequency  calculations,  while  the  upland  extent  of  flood- 
inhibited  rings  may  yield  magnitude  information. 
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The  year-to-year  sequence  of  ring  widths  can  be 
expressed  as  a  time  series  for  statistical  analysis. 
Standardization  techniques  enable  combination  of  individual 
time  series  into  an  averaged  chronology  representative  of  a 
number  of  trees.  Then,  statistical  techniques  such  as 
correlation,  principal  components,  and  regression  can  be 
applied  to  determine  possible  causal  relationships  between 
tree- ring  growth  and  environmental  variables  (Fritts,  1976) . 

The  tree- ring  time  series  can  be  used  to  reconstruct  the 
discharge  time  series,  and  patterns  or  trends  can  be 
identified.  For  exanple,  time  series  analysis  has  shown  that 
California  has  experienced  fewer  extreme  droughts  during  the 
twentieth  century  than  in  previous  centuries  (Hughes  and 
Brown,  1992) .  Tree- ring  chronologies  developed  from  cores  of 
big- cone  spruce,  used  to  reconstruct  precipitation 
variability,  show  major  fluctuations  of  variability  over  the 
past  few  thousand  years.  The  chronologies  showed  the  modern 
period  of  1920-1965  to  be  relatively  dry  with  low 
variability,  with  signs  that  conditions  may  be  becoming  more 
variable  again  (Michaelsen,  et  al.,  1987). 

Despite  the  fact  that  dendrochronology  has  impressive 
dating  accuracy  compared  to  other  methods  of  paleo- study, 
there  are  certain  limitations  that  must  be  considered  when 
interpreting  results. 


Under  maximum  stress  conditions,  tree  growth  for  a 
particular  year  largely  depends  on  that  year's  moisture 
input.  More  commonly,  a  combination  of  soil  moisture  and 
biological  "carry-overs"  result  in  a  lag  effect  (auto¬ 
correlation)  ,  with  a  particular  year's  growth  being 
influenced  by  moisture  from  the  previous  one  or  two  years 


(Smith  and  Stockton,  1981) .  Time  series  analyses  show  that  a 
one  to  two  year  lag  correlates  best  with  hydrologic  discharge 
records  (wendland  and  Watson- Stegner,  1983) .  Autocorrelation 
can  create  uncertainty  in  the  flood  information  stored  in 
tree  rings.  Although  one  can  account  for  a  growth  lag  in  the 
current  data,  it  cannot  be  assumed  that  this  lag  will  apply 
in  the  same  manner  throughout  the  reconstructed  hydrology. 
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CHAPTER  2 
STUDY  ARBA 


2. A  Basin  description 

The  trees  studied  are  located  on  floodplain  and  bluffs 
east  of  the  Big  Muddy  River  where  it  flows  through  the 
Shawnee  National  Forest  (Figure  5) .  Draining  an  area  of  5618 
km2,  the  Big  Mhddy  originates  in  Jefferson  County,  Illinois, 
and  flows  southwestward  for  248  km  to  join  the  Mississippi 
River  near  Grand  Tower  (Illinois  Environmental  Protection 
Agency,  1988) .  The  study  site  for  this  research  is  located 
about  12  km  upstream  of  the  Big  Muddy-Mississippi  confluence. 
The  lower  30  km  of  the  Big  Muddy,  which  includes  the  study 
site,  hugs  the  east  wall  of  a  wide  valley  until  it  sharply 
turns  west  toward  the  Mississippi  River. 

Because  of  its  gentle  channel  slope  and  sluggish 
currents,  the  Big  Muddy  rises  slowly  and  stays  in  flood  for 
relatively  long  periods  (US  Arnty  Corps  of  Engineers,  1968). 
Flooding  on  the  main  stem  of  the  river  is  currently  monitored 
by  nine  stream  gages.  The  nearest  stream  gage  to  the  study 
site  is  located  10  km  upstream,  at  Murphysboro. 

There  are  twelve  reservoirs  of  varied  capacity  in  the 
Big  Muddy  River  drainage  basin,  the  largest  of  which  is  Rend 
Lake.  The  Rend  Lake  reservoir  and  dam,  completed  in  October 
1970,  lies  at  the  northern  end  of  the  watershed.  It  was 
designed  to  reduce  downstream  flooding  by  35%  (U.S.  Army 
Engineer  District,  1970) .  Since  construction  of  the  lake, 
mean  discharge  at  Murphysboro  has  been  reduced  by  5%  (USGS, 
1991) .  The  remainder  of  the  reservoirs  are  intended  for 
recreational  and  water  supply  use,  and  provide  only  limited 
flood  control  (USACOE,  1968;  USAED,  1970). 


Murphy aboro 


Big  Muddy  River  drainage  basin  in  southern 


2.B  Geology  and  Geomorphology 


The  Big  Muddy  basin  lies  between  a  large  sandstone  and 
limestone  cuesta  to  the  east  and  the  Ozark  Plateau  to  the 
west  (Harris,  et  al.,  1977).  The  Shawnee  Hills  region  is  a 
complex,  dissected  upland  located  south  of  past  glaciation 
limits  (USACOE,  1968) .  Glacial  outwash  on  the  Mississippi 
River  channel  created  a  broad,  thick  floodplain  which 
encompasses  the  lower  70  km  of  the  Big  Muddy  River  (USAED, 
1970) .  Bottomland  soils  are  mainly  silt,  clay,  and  loess, 
with  poor  drainage  (US  Department  of  Agriculture,  1975) . 

During  flooding,  sediment  is  carried  out  of  river 
channels  and  onto  the  floodplain.  Coarse  material  is  often 
deposited  near  the  bank,  building  natural  levees  (Harris,  et 
al.,  1977).  The  active  Big  Muddy  floodplain  is  not  very  flat 
as  it  contains  overflow  scour  channels  and  raised  bars. 

Older,  higher  levels  are  rarely  flooded,  while  lower 
surfaces  are  more  frequently  scoured.  Terraces,  at  the 
highest  level  above  the  flat  floodplain,  are  themselves 
remnants  of  an  abandoned  floodplain  (Harris,  et  al.,  1977). 

2.C  Climate  and  vegetation 

The  climate  of  the  lower  Big  Muddy  River  basin  is 
characterized  by  hot  and  humid  summers,  and  relatively  mild 
winters  (Table  1).  Mean  annual  precipitation  is  111.2  cm, 
with  a  peak  during  the  spring.  The  growing  season  for  the 
study  section  of  the  river  lasts  about  200  days,  from  mid- 
April  to  late  October  (USACOE,  1968) .  Summer  precipitation 
is  generally  discontinuous  in  the  form  of  rainshowers,  while 
winter  precipitation  is  more  uniformly  distributed  across  the 
region,  winter  precipitation,  infrequently  snow,  results  in 
a  good  accumulation  of  soil  moisture  by  spring  and  minimizes 
drought  during  slimmer  on  most  soils  in  most  years  (Herman, 
1979) . 
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26 

24 
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Mean 

P  (mm) 

94 

102 

76 

64 

99 

97 

Table  l.  Climate  data  for  Carbondale,  Illinois.  Data  listed 
are:  1961-1991  mean  monthly  values  of  maximum  temperature, 
minimum  temperature,  and  precipitation  (from  Illinois  State 
Water  Survey) . 
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Vegetation  in  the  Shawnee  Hills  is  topographically 
graded  from  bottomland  to  ridge  crest.  Hillslopes  and  ridges 
are  dominated  by  short  grasses  and  several  species  of  oak, 
while  bottomlands  have  a  great  variety  of  flora  to  include 
cottonwood,  maple,  elm,  and  oak  (USACOE,  1968) . 

2.D  Field  site  description 

The  field  site  for  this  research  is  on  floodplain  and 
uplands  east  of  the  Dillon  Bend  meander  of  the  Big  Muddy 
River,  12  km  upstream  from  its  confluence  with  the 
Mississippi  (Figure  6) .  Dillon  Bend  is  bordered  by  100  m 
bluffs  to  the  east  and  abandoned  Mississippi  River  floodplain 
to  the  west.  The  location  of  the  study  site  is  remote, 
accessible  by  unimproved  road  at  the  base  of  the  bluffs  along 
which  there  are  no  habitations  or  permanent  structures.  Land 
surrounding  Dillon  Bend  is  devoted  mainly  to  national  forest, 
although  there  are  a  few  small  agricultural  plots  nearby.  A 
levee  parallels  the  river  about  2  Jan  west  of  Dillon  Bend 
(Figure  7).  The  levee,  part  of  the  Grand  Tower  Drainage  and 
Levee  District,  was  completed  in  1959  to  protect  residents 
and  agricultural  land  from  high  flows  on  the  Big  Muddy  and 
Mississippi  (Dobney,  1975) . 

The  landscape  at  the  Dillon  Bend  study  site  lias  a 
definite  lower  floodplain  level  interrupted  by  higher  banks 
and  natural  levees.  The  banks  and  levees  vary  in  width; 

5-10  m  in  some  spots  while  spanning  broad  fields  in  other 
locations.  If  the  banks  are  remnant  terraces  or  ancient 
river  sediments,  some  soil  development  would  most  likely  have 
occurred  (Johnson,  personal  communication,  1992;  Martin, 

1992)  .  A  natural  cutlaank  exposure  cut  into  the  fill  of  a 
terrace  was  examined  for  soil  development.  The  top  of  the 
cutbank  is  about  4  m  above  the  bankfull  level  of  the  Big 
Muddy  River.  The  t>ank  contains  unstratified  fill  with  signs 
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Figure  6.  Dillon  Bend  section  of  the  Big  Muddy  River  (after 
USGS,  1990) . 
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of  leaching,  but  with  no  visible  sign  of  a  clay  (B)  horizon. 
This  suggests  that  the  highest  surface  on  the  eastern  side  of 
Dillon  Bend  is  relatively  young  and  aggrading,  and  most 
likely  a  historically  active  floodplain. 

2.1  Known  flood  history 

The  floodplain  along  the  Dillon  Bend  section  of  the  Big 
Muddy  River  is  prone  to  flooding  under  two  separate 
scenarios.  The  first  situation  involves  flows  resulting  from 
runoff  within  the  Big  Muddy  drainage  basin,  and  is 
characterized  by  high  gaged  discharge  at  the  Murphysboro 
stream  gage.  Secondly,  high  flows  on  the  Mississippi  River 
may  result  in  backwater  flooding  of  the  lower  60  to  80  km  of 
the  Big  Muddy,  which  includes  Dillon  Bend  (USAED,  1970) .  The 
nearest  streamgage  to  the  Big  Muddy/Mississippi  River 
confluence  is  at  Thebes,  Illinois,  51  km  downstream  from 
Dillon  Bend. 

The  occurrence  and  degree  of  flood  events  is  best 
measured  with  values  of  peak  stage,  or  the  maximum  height  of 
floodwaters  at  a  gaging  station.  The  known  flood  history, 
with  related  stages,  for  the  lower  Big  Muddy  River  from  gaged 
records  and  historical  data  can  be  summarized  as  follows: 

2.E.1 _ Flooding  from  high  Lines,  of  the  .Big. Muddy-  Slyer 

There  have  been  numerous  historical  reports  of  damaging 
floods  in  the  Big  Muddy  River  basin  since  the  early  1900s. 
"Newsworthy"  flood  events  were  all  associated  with  peak 
stages  of  at  least  10.0  m  at  the  Murphysboro  gage  (USACOE, 
1968;  USAED,  1970).  The  Murphysboro  gage  has  been  recording 
river  stage  since  1916,  but  data  is  fragmentary  prior  to 
1931.  Significant  flooding  was  reported  by  area  residents  in 
March  1913  and  August  1915,  but  nothing  is  known  about  the 
magnitudes  of  these  events. 
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From  1916  to  present,  the  reported  stage  at  Murphysboro 
has  reached  10. 0  m  or  greater  11  times  (Table  2).  The  flood 
of  record  occurred  on  May  12,  1961,  when  the  stage  peaked  at 
11.6  m.  A  major  storm  over  the  basin  had  produced  widespread 
rain  and  the  record  discharge  for  all  gaging  stations  on  the 
main  stem  (USACOE,  1968).  Relative  magnitudes  of  other  gaged 
flood  events  are  compared  with  values  of  peak  stage  in 
Table  2. 

2.E.2 _ Backwater  flooding  from  the  Mississippi  .River 

Comparison  of  gaged  records  between  Murphysboro  and 
Thebes  shows  that  high  stages  usually  do  not  occur 
simultaneously  at  the  two  locations.  Thebes  discharge 
includes  inputs  from  the  entire  Upper  Mississippi  River 
watershed  (1,847,188  lan2),  while  flow  at  Murphysboro  responds 
to  runoff  from  its  own,  much  smaller  watershed  (5,618  km2)  . 
Major  floods  on  the  Upper  Mississippi  River  usually  occur 
between  March  and  June,  from  rapid  melting  of  a  large 
accumulation  of  snow  (USAED,  1970)  .  Considering  that  weather 
systems  and  climatic  conditions  that  contribute  to  runoff  and 
snowmelt  vary  dramatically  over  such  a  large  area,  and  that 
there  is  a  significant  lag  time  for  Thebes  discharge  to 
respond  to  runoff  inputs,  it  is  not  surprising  that  peak 
stages  between  Thebes  and  Murphysboro  are  not  coincident. 

The  streamgage  at  Thebes  has  been  recording  discharge 
data  since  1934,  but  stage  measurements  from  an  auxiliary 
gage  13  km  upstream  date  back  to  1896.  The  only  known  major 
flood  outside  the  period  of  record  occurred  in  July,  1844, 
with  an  estimated  stage  of  13.0  m,  computed  by  the  U.S.  Army 
Corps  of  Engineers  using  floodmarks  (USGS,  1991) .  The 
Mississippi  River  floods  of  1943,  1944,  and  1947  are 
documented  to  have  caused  backwater  flooding  on  the  lower 
section  of  the  Big  Muddy  River,  contributing  to  the  high 
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Flood- Stage  Data  for  Murphysboro,  Illinois 


ElQQfl _ 

May  1933 
March  1935 
January  1937 
August  1946 
January  1949 
January  1950 
April  1957 
May  1961 
March  1979 
May  1983 
April  1985 


Peak  Stage  at  Murphysboro  (m) 

10.3 

10.4 
10.3 

10.1 
11.0 
10.8 

10.3 

11.6 

10.4 

11.2 

10.2 


Table  2.  Stage  data  for  major  floods  (greater  than  10. 0  m) 
in  the  Big  Muddy  River  basin,  recorded  at  the  Murphysboro 
gaging  station  (from  USACOE,  1968;  united  States  Geological 
Survey,  1991) . 
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stages  recorded  at  Murphysboro  for  those  years  (USACOE, 

1968) .  A  farmer  who  currently  manages  fields  along  the 
eastern  bank  of  the  Big  Muddy  south  of  Dillon  Bend,  confirmed 
that  the  area  is  highly  prone  to  backflooding  from  the 
Mississippi  River.  The  backwaters  are  said  to  overflow  the 
banks  near  Dillon  Bend  when  the  stage  at  the  Thebes 
streamgage  measures  about  9.7  m  (Dennis  Ballance,  pers. 
comm.,  1992). 

Since  1932,  the  stage  of  the  Mississippi  River  at  Thebes 
has  exceeded  9.7  m  47  times  in  45  years.  Peak  stage  on 
record  at  Thebes  occurred  with  the  April  1973  flood,  when  the 
floodwaters  reached  13.2  m.  The  rare  event  occurred  as  the 
Mississippi,  Missouri,  and  Illinois  Rivers  crested 
simultaneously  (Dobney,  1978) .  This  and  other  floods  are 
compared  in  Table  3. 
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Flood- Stage  Data  for  Thebes,  Illinois 


Flood 

Peak  Sfcaae  (m) 

Flood  Peak 

Stage  (m) 

July  1844 

13.0* 

July  1951 

12.1 

June  1903 

11.1 

May  1952 

11.1 

May  1904 

10.4 

July  1958 

10.3 

June  1908 

10.4 

April  1960 

11.3 

July  1909 

10.6 

May  1961 

11.8 

April  1912 

10.6 

March  1962 

11.0 

August  1915 

10.4 

April  1965 

10.4 

February  1916 

11.1 

October  1966 

9.9 

June  1916 

11.1 

July  1969 

11.8 

April  1920 

9.8 

May  1970 

11.0 

April  1922 

11.6 

March  1971 

10.0 

October  1926 

10.0 

April  1972 

10.0 

April  1927 

12.2 

April  1973 

13.2 

June  1928 

10.0 

May  1974 

10.0 

May  1929 

11.4 

April  1975 

11.0 

May  1933 

10.5 

March  1978 

11.6 

June  1935 

11.1 

April  1979 

13.1 

June  1942 

10.8 

May  1983 

12.1 

May  1943 

12.2 

April  1984 

11.7 

May  1944 

11.9 

March  1985 

12.3 

April  1945 

11.5 

November  1985 

12.3 

January  1946 

10.2 

October  1986 

12.6 

July  1947 

12.2 

May  1990 

12.1 

May  1950 

9.8 

*  estimated  from  floodmarks 


Tahi a  stage  data  for  major  floods  on  the  Upper 
Mississippi  River,  recorded  at  the  Thebes  gaging  station 
(from  USACOE,  1945;  USQS,  1993) 
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CHAPTER  3 
METHODS 


3. A  Data  acquisition  and  analysis 
3„>AJL  Trea  cores 
Coring  procedures: 

Objectives  at  the  research  site  were  to  sample  trees  of 
a  favorable  species  and  longevity  for  dendroclimatic 
research,  and  to  core  several  levels  above  the  floocfclain  on 
the  hillslope  and  terrace,  u.s.  Forest  Service  personnel, 
Harrisburg  office,  provided  assistance  in  locating  old  tree 
stands.  Oak  and  hackberry  were  sampled  due  to  their 
abundance,  longevity,  and  proven  response  to  hydrologic 
stimuli  (Hosner  and  Boyce,  1962) .  Collected  species  include 
black  oak  (Quercus  velUtina  Lam. ) ,  white  oak  (Ouercus  alba 
L. ) ,  bur  oak  (Ouercus  macrocarpa  Michx. ) ,  chinakpin  oak 
(Ouercus  mmihlenberyi i  Engelm. ) ,  shumard  oak  (Ouercus 
Bhumarfli i  Buckl.),  red  oak  (Ouercus  Rubra  L.),  and  hackberry 
(Celt is  occidentalls  Torr.)  .  Minimum  sample  size  at  each 
elevational  level  was  five  trees  for  chronology  development. 
Sampling  distribution  is  shown  in  Figure  7. 

At  least  two  cores  per  tree,  separated  by  90°,  were 
extracted  to  account  for  variable  and  irregular  growth 
patterns  within  the  stem  (Stokes  and  Smiley,  1968) .  Parts  of 
the  stem  where  rings  might  be  distorted  (such  as 
uphill/downhill  sides  and  parts  near  brunches)  were  not 
sampled.  Cores  were  directed  toward  the  pith  to  obtain  a 
full  record  of  the  tree's  longevity. 

Tree  cores  were  taken  with  a  Haglof  4.3  mm  increment 
borer,  which  screws  into  the  tree  and  removes  a  cylindrical 
sample  i.i  which  the  rings  are  preserved.  Cores  were  removed 
from  the  borer  and  placed  in  plastic  straws ,  sealed  and 
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•  White  oak 
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a  Shumard  oak 
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Figure  7  -  Distribution  of  trees  sanpled  along  the  Dillon 
Bend  reach  of  the  Big  Muddy  River.  Horizontal  distances  are 
not  to  scale.  Spot  elevations  represent  surveyed  ground 
level  in  meters  above  sea  level  (MSL) . 
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numbered.  Notes  were  taken  on  location,  species,  elevation 
above  the  river,  core  orientation,  and  physical  condition  of 
the  tree  (Stokes  and  Smiley,  1968) .  Each  sampled  tree  was 
labeled  with  a  metal  tag. 

Ring  measurement: 

Cores  were  allowed  to  dry  for  a  few  days  in  the  straws 
to  prevent  warping.  They  wer'*  then  mounted  on  pine  strips 
and  sanded  with  medium  grade  sandpaper  to  give  a  flat, 
measurable  surface.  Sanded  surfaces  were  shaved  using  a 
razorblade  to  expose  cells,  and  given  a  final  sanding  with 
400  grade  sandpaper.  After  the  rings  were  counted  and  marked 
by  decades,  the  width  of  each  annual  ring  was  measured  and 
recorded  using  a  stage  microscope.  For  consistent 
methodology,  all  measurements  were  made  perpendicularly 
across  rings  from  latewood  to  earlywood  boundaries  so  that 
the  relative  ring  widths  could  be  compared  (Stokes  and 
Smiley,  1968) . 

Flood  ring  identification: 

Low  level  (floodplain)  and  mid- level  cores  were 
inspected  under  a  microscope  for  wood  abnormalities  related 
to  flooding.  Flood  rings  were  identified  and  annotated  for 
growth  year  and  intra-ring  position  of  the  abnormal  growth 
(i.e.  mid- latewood) .  Sing  position  can  be  used  to  estimate 
the  relative  time  of  flooding  during  the  growing  season. 

After  inspection  of  cores  was  completed,  flood  ring  years 
were  marked  on  a  time  chart  with  separate  time  lines  for  each 
tree,  so  that  ccercnonalities  could  be  identified. 

3. A. 2  Hydrologic  data 

Stream  gage  records  provide  mean  daily  discharge  data 
(cubic  meters  per  second)  and  extreme  peak  discharges  (m3/s) 
for  the  duration  of  a  station's  operation.  Mean  monthly 
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discharge  data  from  two  gages  were  collected  with  assistance 
from  personnel  at  United  States  Geological  Survey,  Urbana 
office.  Gages  at  Thebes  and  Murphysboro  were  selected  for 
proximity  to  the  study  site  and  length  of  record.  The  Thebes 
gage  monitors  the  Mississippi  River  51  km  downstream  from  the 
Mississippi/Big  Muddy  confluence,  recording  from  1934  to 
present.  The  Murphysboro  gage  is  located  on  the  Big  Muddy 
River  10  km  upstream  from  Dillon  Bend,  and  has  been  operating 
since  1916  (fragmentary  prior  to  1931) (USGS,  1991) . 

3. A. 3 _ Climate  data 

The  National  Weather  Service  maintains  daily  records  of 
temperature  and  precipitation  for  some  200  locations  in 
Illinois.  For  this  study,  average  maximum  temperature  and 
precipitation  totals  for  Carbondale  and  Anna  were  collected 
for  each  month  of  each  year  since  1901.  The  data  can  be  used 
to  determine  the  effect  of  temperature  and  precipitation  on 
tree  growth  during  the  growing  season,  and  during  the  spring 
months  when  growth  begins.  Precipitation  during  the  fall  and 
winter  months  prior  to  the  next  growing  season  can  be  used  to 
study  the  influence  of  soil  moisture  storage. 

3.B  Model  development 

3.B.1  Craasdatiag 

Before  continuing  with  model  development,  it  is 
necessary  to  accurately  date  the  rings  within  the  tree  cores. 
This  includes  accounting  for  missing  rings  at  the  ends  of  or 
within  the  cores,  and  multiple  or  false  rings.  Since  trees 
record  the  environment  in  a  pattern  of  wide  and  narrow  rings, 
it  is  possible  to  match  patterns,  or  crossdate  trees  (Fritts, 
1976) .  The  calendar  year  of  the  outermost  ring,  known  when 
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bark  is  present,  serves  as  the  starting  point  for 
crossdating. 

The  procedure  of  crossdating  uses  the  COFECHA  computer 
program  (see  Appendix  B) ,  which  aids  in  identifying  errors 
and  acts  as  an  independent  tool  to  confirm  the  accuracy  of 
dating  and  measurement  (Holmes,  et  al.,  1986).  The  dated  and 
measured  ring  series  are  filtered  by  fitting  a  20 -year  cubic 
spline,  and  then  dividing  the  series  values  by  the 
corresponding  spline  curve  values  to  remove  low-frequency 
variance.  The  high  frequency  residual  is  subjected  to  a  log 
transformation  to  equalize  proportionally  the  variability 
among  small  and  large  rings  (Holmes,  et  al.,  1986).  The  mean 
values  of  the  combined  series  comprise  the  master  dating 
chronology.  After  an  initial  COFECHA  run  using  the  40  oak 
ring  series,  the  10  series  with  the  best  correlation  to  the 
master  are  selected  to  compose  a  more  reliable  master  dating 
chronology.  Each  COFECHA  output  identifies  measuring  or 
counting  errors  by  flagging  40  year  ring  segments  and 
individual  outliers  within  each  core  that  do  not  agree  with 
the  master.  After  using  the  output  to  account  for  possible 
problems  in  master  cores,  each  of  the  remaining  cores  are 
added  one  by  one  to  the  COFECHA  program  and  checked  against 
the  master  chronology.  As  measuring  or  counting  errors  are 
remedied,  each  core  is  included  in  the  final  dating 
chronology . 

3.B.2 _ Standardization 

Ring-widths  generally  decrease  as  trees  age,  due  to  the 
tree's  ever-increasing  diameter  (Fritts,  1976).  To  compare 
trees  with  different  growth  rates  and  ages,  growth  trend 
analysis  procedures  must  be  used  to  produce  annual 
standardized,  dimensionless  indices  for  all  rings.  The 
computer  program  ARSTAN  (see  Appendix  b)  fits  either  negative 
exponential  curves  or  stiff  cubic  splines  to  tree- ring  series 
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(Holmes,  et  al.,  1986).  Trend  curves  differ  in  form 
depending  upon  site- specif ic  growth  variations  and  distance 
from  the  pith  of  the  tree,  visual  inspection  of  raw  ring- 
widths  is  used  to  estimate  whether  the  trend  of  growth  is  a 
clear  departure  from  negative  exponential.  ARSTAN  prompts  as 
to  whether  the  operator  believes  there  are  any  cores  with 
clear  departures.  Specified  cores  are  fitted  with  cubic 
splines  to  model  growth  trend.  Figures  8-10  show  raw  ring- 
widths  plotted  against  their  computer- fitted  growth  curves. 
Note  that  cores  DB37,  DB38,  and  DB43  did  not  have  negative 
exponential  growth  trends,  requiring  cubic  spline  treatment. 

An  annual  index  is  computed  by  dividing  the  measured 
ring  width  by  the  value  of  the  growth  curve  for  that 
particular  year.  To  account  for  the  lag  influence  of 
previous  years'  growth,  ARSTAN  uses  autoregressive  modeling 
to  enhance  the  common  signal  (Holmes,  et  al.,  1986).  The 
chronology  of  a  tree  core  is  then  expressed  as  a  series  of 
indices  that  are  stationary  in  both  mean  and  variance 
(Stockton,  1975) . 

Master  chronologies  are  constructed  for  the  floodplain, 
mid- slope,  and  high  level  trees  by  combining  the  indices  from 
individual  cores  from  each  level,  and  computing  an  average 
growth  index  for  each  year  in  the  combined  chronology.  Each 
of  the  three  levels  are  then  represented  by  a  single  master 
chronology. 

3.B.3 _ Correlation  analysis 

Paleoflood  analysis  requires  testing  of  the  relationship 
between  tree- ring  growth,  climate,  and  discharge.  Prior  to  a 
regression  analysis,  it  is  necessary  to  measure  correlation 
between  these  variables,  so  that  significant  months  or 
seasons  (groupings  of  months)  can  be  identified.  Pearson 
correlation  coefficients  with  a  95%  level  of  confidence  are 
identified  between  environmental  variables  and  the  following 
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^  Raw  ring-widths  from  low-level  (floodplain)  trees  fitted  with  growth  trend 
using  the  computer  program  ARSTAN.  Core  species  are  bur  oak,  except  for  DB44 


Dillon  Bend  Mid-Level  Site 
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Figure  9 .  Raw  ring-widths  from  mid- level  trees  fitted  with  growth  trend  curves  using 
the  computer  program  ARSTAN.  Core  species  for  DB15  through  DB20  are  white  oak.  DB37 
and  DB38  are  black  oak.  DB25  is  bur  oak. 
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Figure  IQ.  Raw  ring-widths  from  high-level  trees  fitted  with  growth  trend  curves  using 
the  computer  program  ARSTAN.  Core  species  are  black  oak:  DB39,  DB40,  DB49,  DB50,  DB53, 
DB54;  bur  oak:  DB14,  DB14;  chinkapin  oak:  DB51,  DB52;  white  oak:  DB21,  DB22,  DB41, 

DB42;  shumard  oak:  DB5 5,  DB5 6,  DB57 . 


growth  variables:  floodplain  trees,  mid- slope  trees,  high 
level  trees,  and  ratios  of  floodplain/high  trees. 
Coefficients  are  calculated  for  all  12  months  of  the  growth 
year,  and  May  through  December  of  the  previous  year  to 
account  for  lag  effects.  With  tree  growth  as  the  dependent 
variable,  independent  variables  were  selected  for  the 
following  reasons: 

a.  Carbondale  monthly  precipitation  and  temperature 

Determines  the  degree  of  a  linear  relationship  between 
tree- growth  and  local  precipitation  and  temperature 
recorded  at  the  nearest  climate  station  within  the  Big 
Muddy  watershed. 

b.  Anna  monthly  precipitation  and  temperature 

Determines  the  degree  of  a  linear  relationship  between 
tree-growth  and  local  precipitation  and  temperature 
recorded  at  the  nearest  weather  station  in  an  adjacent 
watershed. 

c.  Regional  precipitation  and  temperature  for  Illinois 
Climate  Division  8. 

Determines  the  degree  of  a  linear  relationship  between 
tree-growth  and  precipitation  and  temperature  from  the 
lower  section  of  the  Upper  Mississippi  watershed  above 
the  Mississippi-Ohio  River  confluence. 

d.  Illinois  regional  Palmer  Drought  Severity  Index  (PDSI) . 

Determines  the  degree  of  a  linear  relationship  between 
tree  growth  and  drought  as  defined  by  PDSI  (see 
Appendix  A) . 
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e.  Illinois  regional  Palmer  Hydrologic  Drought  Index  (PHDI) . 

Determines  the  degree  of  a  linear  relationship  between 
tree  growth  and  soil  moisture  as  defined  by  PHDI  (see 
Appendix  A)  . 

f.  Murphysboro  mean  monthly  discharge. 

Determines  the  degree  of  a  linear  relationship  between 
tree -growth  and  flow  from  the  Big  Muddy  River  watershed. 

g.  Thebes  mean  monthly  discharge. 

Determines  the  degree  of  a  linear  relationship  between 
tree  growth  and  flow  from  the  Mississippi  River 
watershed  north  of  Thebes  (known  to  backflood  the  Big 
Muddy  channel) . 

After  correlations  are  calculated  for  each  of  these 
monthly  variables,  relationships  which  are  significant  at  the 
0.05  level  are  identified.  Consecutive  significant  months 
are  grouped  into  "seasons",  and  correlations  between  these 
aggregate  variables  and  tree-ring  growth  are  calculated. 
Further  statistical  analysis  uses  the  seasonal  variables  with 
the  highest  correlations,  provided  that  they  are  reasonable 
environmental  parameters. 

3 -B.4..  .Comparisons 

To  investigate  the  effects  of  flooding  on  floodplain 
trees  compared  to  upland  trees,  ratios  of  floodfc>lain/upland 
standardized  ring-width  indices  are  computed  for  each  year  of 
the  tree-ring  chronology  (see  l.D.3).  The  ratios  are  first 
plotted  against  normals  ed  seasonal  discharge  data  to 
identify  how  well  trends  in  tree- ring  ratios  conform  to 
trends  in  discharge.  Then,  the  best - correlated  discharge 
variable  (determined  from  correlation  analysis)  is  regressed 
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on  the  ratio  variable  to  develop  a  predictive  model  for  flood 
reconstruction . 

3.B..5 _ Bivariate  linear  regression 

Since  the  nature  of  this  research  is  flood 
reconstruction,  it  is  necessary  to  analyze  the  relationship 
between  discharge  and  tree- ring  growth.  This  is  accomplished 
by  regressing  discharge  on  the  standardized  chronologies, 
where  discharge  variables  are  the  seasonal  totals  selected 
from  correlation  analysis.  One  half  of  the  discharge  data 
set  is  used  for  model  calibration,  and  the  other  saved  for 
verification  of  the  derived  model.  The  SAS  statistical 
software  package  is  used  to  compute  model  parameters  and 
analysis  of  variance  (ANOVA)  test  statistics.  The  regression 
equation  is  verified  against  the  remaining  half  data  set 
using  the  verify  software  package  (see  Appendix  B)  (Holmes, 
et  al.,  1986).  Verification  results  are  tested  at  the  95% 
confidence  level  for  correlation,  reduction  of  error,  and  T- 
value. 
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CHAPTER  4 
RESULTS 


4. A  Chronology  development 

Two  of  the  most  important  procedures  in  chronology- 
development  are  crossdating  and  standardization .  Crossdating 
assures  proper  placement  in  time  of  each  annual  ring,  and 
standardization  corrects  for  systematic  growth  changes 
resulting  from  tree  aging  (Fritts,  1976) .  Separate 
chronologies  were  developed  for  the  floodplain,  mid- slope  and 
high  level  collection  sites  at  Dillon  Bend.  Specifics  on 
chronology  development  are  explained  in  this  section. 

1..A.1 _ Crossdatinq  and  verification 

The  final  floodplain,  mid- slope  and  high  level  data  sets 
consist  of  12,  9,  and  17  cores,  respectively.  The  computer 
program  COFECHA  (Holmes,  et  al.,  1986)  was  used  to  identify 
and  correct  possible  dating  or  measurement  errors  in  the  38 
cores  (see  Appendix  B) .  A  total  of  4  cores  were  emitted  from 
the  data  set  because  of  poor  correlation  with  the  master 
series  for  their  respective  site.  Three  of  the  4  abandoned 
cores  were  from  the  floodplain  site,  which  has  the  most 
variable  growth  regime  due  to  differential  drainage  from 
uneven  topography. 

The  minimum  acceptable  correlation  coefficient  in 
crossdating  was  0.3665,  which  represents  a  99%  level  of 
confidence  when  testing  40  year  segments  of  the  cores.  A 
segment  length  of  40  years  provided  sufficient  degrees  of 
freedom  to  avoid  chance  correlations,  while  being  short 
enough  to  help  locate  dating  problems  (Holmes,  et  al.,  1986). 
Descriptive  statistics  on  the  final  chronology  for  each  site 
are  listed  in  Tables  4-6. 
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Number  of  dated  series:  12 

Master  series  1905-1991:  87  yrs 

Average  mean  sensitivity:  .221 

Series  intercorrelation:  .538 


■  Mil  '  >1  III  WM 

_ CQEr.— W/-J 

DB01 

1916*1991  76 

0 

.588 

DB02 

1921*1991  71 

0 

.524 

DB05 

1915*1991  77 

1 

.539 

DB06 

1915-1991  77 

0 

.552 

DB07 

1923-1979  57 

0 

.481 

DB08 

1931-1979  49 

0 

.473 

DB09 

1905-1991  87 

0 

.605 

DB10 

1905-1939  35 

0 

.585 

DB11 

1939-1991  53 

0 

.466 

DB12 

1929-1991  63 

0 

.494 

DB43 

1920-1991  72 

0 

.530 

DB44 

1952-1991  40 

0 

Mean:  .538 

Tahi«a  a.  Descriptive  statistics  on  floodplain  tree  cores 
after  completion  of  crossdating  procedures.  #  Flags  refers  to 
the  number  of  40  year  segments  which  had  less  them  .3665 
correlation  with  the  master  floodplain  chronology.  The 
overall  series  correlation  with  the  floodplain  master  is 
listed  in  the  final  column. 


DILLON  BEND  MID- SLOPE  TREES 


Number  of  dated  series:  9 

Master  series  1892-1991:  100  yrs 

Average  mean  sensitivity:  .239 

Series  intercorrelation:  .664 


DB15 

1924-1991 

68 

0 

.597 

DB16 

1924-1991 

68 

0 

.630 

DB17 

1892-1991 

100 

0 

.711 

DB18 

1901-1991 

91 

0 

.721 

DB19 

1906-1991 

86 

0 

.720 

DB20 

1912-1991 

80 

0 

.737 

DB25 

1917-1991 

75 

1 

.530 

DB37 

1950-1991 

42 

0 

.663 

DB38 

1925-1991 

67 

1 

.617 

Mean:  .664 


Table  5 .  Descriptive  statistics  on  mid- slope  tree  cores 
after  completion  of  crossdating  procedures.  #  Flags  refers  to 
the  number  of  40  year  segments  which  had  less  than  .3665 
correlation  with  the  master  mid- slope  chronology.  The 
overall  series  correlation  with  the  mid- slope  master  is 
listed  in  the  final  column. 
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Number  of  dated  series:  17 

Master  series  1905-1991:  123  yrs 

Average  mean  sensitivity:  .220 

Series  intercorrelation:  .643 


H  ill 

il  l' 

■l¥4  'll  — 

'JM  ».  wThi 

—  1  m 

DB13 

1912-1979 

68 

0 

.552 

DB14 

1912-1979 

68 

0 

.623 

DB21 

1882-1979 

98 

0 

.510 

DB22 

1869-1979 

111 

1 

.514 

DB39 

1921-1991 

71 

0 

.643 

DB40 

1935-1989 

55 

0 

.547 

DB41 

1898-1991 

94 

0 

.610 

DB42 

1885-1939 

107 

0 

.666 

DB49 

1922-1991 

70 

0 

.676 

DB50 

1915-1991 

77 

0 

.762 

DB51 

1912-1991 

60 

0 

.694 

DB52 

1909-1959 

51 

0 

.607 

DB53 

1904-1991 

88 

0 

.685 

DB54 

1912-1991 

80 

0 

.720 

DB55 

1914-1991 

78 

0 

.631 

DB56 

1905-1991 

87 

0 

.799 

DB57 

1945-1991 

47 

0 

.718 

Mean:  .643 


Table  6.  Descriptive  statistics  on  high-level  tree  cores 
after  completion  of  crossdating  procedures.  #  Flags  refers  to 
the  number  of  40  year  segments  which  had  less  than  0.3665 
correlation  with  the  master  high-level  chronology.  The 
overall  series  correlation  with  the  high-level  master  is 
listed  in  the  final  column. 


The  mean  correlation  coefficients  for  the  floodplain 
(0.538),  mid-slope  (0.664),  and  high  sites  (0.643)  are  very 
good  for  a  deciduous  species  such  as  oak  (Estes,  1970;  Henri 
Grissino- Mayer,  pers.  comm.,  1992).  It  is  not  surprising 
that  the  mid- slope  and  high  sites  have  higher  mean 
correlations  than  the  floodplain  site,  since  their  climate 
sensitive  landscape  position  favors  a  common  response  among 
trees. 

Another  important  computed  statistic  is  mean 
sensitivity,  which  is  a  measure  of  change  in  ring  width  from 
one  year  to  the  next.  The  average  mean  sensitivity  in  many 
dendrochronologic  studies  is  0.35,  with  a  99%  confidence 
interval  of  0.13  to  0.57  (Fritts,  1991).  The  Dillon  Bend 
chronologies  have  mean  sensitivities  ranging  from  0.220  to 
0.239,  well  within  common  limits.  The  value  of  0.221  for  the 
floodplain  site  was  surprisingly  high  when  considering  its 
poor  drainage  and  abundant  moisture  supply.  This  unexpected 
result  may  reflect  the  fact  that  the  floodplain  is  uneven, 
with  banks  and  scour  holes  which  vary  the  relief,  drainage 
pattern,  and  effects  of  floodwaters. 

Before  finalizing  a  chronology,  it  is  standard  practice 
to  verify  core  measurement  and  dating  by  having  a  second 
party  remeasure  10%  of  the  cores  which  conpose  a  chronology 
(Fritts,  1976) .  Henri  Grissino- Mayer  from  the  Laboratory  of 
Tree-Ring  Research  at  the  University  of  Arizona  randomly 
selected  4  of  the  38  oak  cores  and  remeasured  than  using 
equipment  at  the  University  of  Arizona,  Tucson.  He  verified 
his  measurements  against  my  own  using  the  MEASURE  computer 
program  (see  Appendix  B) .  Inspection  of  descriptive 
statistics  (average  ring  width,  standard  deviation,  etc.)  for 
the  original  and  remeasured  cores  suggests  an  accurate 
measurement,  which  was  then  confirmed  by  a  Matched  Pairs 
difference  of  means  test  at  the  95%  level  of  confidence  (see 
Tables  7-10) . 
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VERIFICATION  RESULTS 


Core  selected:  DB01 


Descriptive  statistics 


Original  file _ Verification  file 


Average  width  (mm) 

4.027 

4.028 

Median  (mm) 

3.990 

3.975 

Variance 

1.956 

1.992 

Std.  deviation  (nsn) 

1.399 

1.411 

Coef.  variation 

0.347 

0.350 

Mean  sensitivity 

0.205 

0.192 

Autocorrelation 

0.746 

0.766 

Regression  results 

Dependent  variable  (y)  original  DB01  file 
Independent  variable  (x)  verification  file 
Slope  of  regression  0.984 

Intercept  of  regression  (mm)  0.064 

Coefficient  of  variation  4.165 

Correlation  coefficient  0.993 


Verification  results:  1916-1991 


Test  for  Matched  Pairs 


Sum  of  differences 

9.200 

Total  square  differences 

2.121 

Avg  absolute  differences 

0.121 

Mean  square  differences 

0.028 

Test  statistic  t 

1.043 

Accepted  at  95%  level  of  confidence 


Table  7 .  Output  statistics  from  the  computer  program 
MEASURE,  which  compares  two  independent  measurements  of  core 
DB01 . 
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VERIFICATION  RESULTS 


Core  selected:  DB08 


Descriptive 

statistics 

Original  file 

Verification 

Average  width  (mm) 

4.321 

4.321 

Median  (mm) 

4.470 

4.400 

Variance 

2.942 

2.973 

Std.  deviation  (mm) 

1.715 

1.724 

Coef.  variation 

0.397 

0.399 

Mean  sensitivity 

0.243 

0.237 

Autocorrelation 

0.750 

0.750 

Regression  results 

Dependent  variable 

(y)  original  DB08  file 

Independent  variable  (x)  verification  file 
Slope  of  regression  0.991 

Intercept  of  regression  (mm)  0.038 

Coefficient  of  variation  3.426 

Correlation  coefficient  0.996 


Verification  results:  1931-1989 


Test  for  Matched  Pairs 


Sum  of  differences 

6.230 

Total  square  differences 

1.263 

Avg  absolute  differences 

0.106 

Mean  square  differences 

0.021 

Test  statistic  t 

1.060 

Accepted  at  £5$ -level  of.  .confidence 


Table  8.  Output  statistics  from  the  conputer  program 
MEASURE,  which  compares  two  independent  measurements  of  core 
DB06. 
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VERIFICATION  RESULTS 


Core  selected:  DB21 


Descriptive  statistics 


Original  file _ verification  file 


Average  width  (mm) 

1.979 

1.979 

Median  (mm) 

2.040 

2.115 

Variance 

0.540 

0.557 

Std.  deviation  (mm) 

0.735 

0.747 

Coef.  variation 

0.371 

0.377 

Mean  sensitivity 

0.190 

0.208 

Autocorrelation 

0.803 

0.788 

Regression  results 

Dependent  variable  (y)  original  DB21  file 
Independent  variable  (x)  verification  file 
Slope  of  regression  0.961 

Intercept  of  regression  (mm)  0.077 

Coefficient  of  variation  7.937 

Correlation  coefficient  0.977 


Verification  results:  1882-1991 
Test  for  Matched  Pairs 


Sum  of  differences  10.990 
Total  square  differences  2.755 
Avg  absolute  differences  0.100 
Mean  square  differences  0.025 
Test  statistic  t  0.997 


Accepted  at  95%  level  of  confidence 


Table  q.  Output  statistics  from  the  computer  program 
MEASURE,  which  compares  two  independent  measurements  of  core 
DB21 . 
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VERIFICATION  RESULTS 


Core  selected:  DB38 


Descriptive  statistics 
original  file _ veri£ i cation  ..file 


Average  width  (mm)  3.281  3.279 
Median  (mm)  3.350  3.310 
Variance  2.236  2.353 
Std.  deviation  (mm)  1.495  1.534 
Coef.  variation  0.456  0.468 
Mean  sensitivity  0.200  0.215 
Autocorrelation  0.856  0.843 


Regression  results 

Dependent  variable  (y)  original  DB38  file 
Independent  variable  (x)  verification  file 


Slope  of  regression  0.972 
Intercept  of  regression  (mm)  0.095 
Coefficient  of  variation  3.629 
Correlation  coefficient  0.997 


verification  results:  1925-1991 
Test  for  Matched  Pairs 


Sum  of  differences  6.500 
Total  square  differences  1.046 
Avg  absolute  differences  0.097 
Mean  square  differences  0.016 
Test  statistic  t  1.201 


Accepted  at  95%  level  of  confidence 


Table  10.  Output  statistics  from  the  computer  program 
MEASURE,  which  compares  two  independent  measurements  of  core 
DB38 . 
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A  standardized  master  chronology  was  computed  for  the 
floodplain,  mid- slope  and  high  collection  levels  using  the 
ARSTAN  computer  program  (see  Appendix  b) .  Ring  width  data 
for  all  cores  were  entered  into  the  program,  which  applied 
growth  trend  analysis  and  autoregressive  modeling  to  produce 
a  standardized  index  value  for  each  year  at  each  of  the 
collection  levels  (Holmes,  et  al.,  1986).  The  mean  value  of 
the  indices  for  each  chronology  approximated  1.0,  allowing 
data  from  trees  with  different  growth  rates  to  be  combined  or 
compared  for  climatic  variation.  Standardized  indices  for 
each  chronology  are  shown  in  Figure  11  and  numerical  values 
are  listed  in  Appendix  D. 

Verification  of  a  regional  climatic  signal  was 
accomplished  by  comparing  standardized  indices  of  the  Dillon 
Bend  high  elevation  cores  with  those  from  2  nearby  sites,  as 
suggested  by  Fritts  (1991),  Kutzbach  and  Guetter  (1980),  and 
Estes  (1970) .  The  nearest  tree- ring  study  to  the  Dillon  Bend 
area  was  that  of  Estes,  in  the  Pine  Hills  area  east  of  the 
Big  Muddy  River,  about  11  km  south  of  Dillon  Bend.  He 
sampled  shortleaf  pine  for  correlation  with  precipitation 
data.  Biasing  and  Duvick  (1984)  collected  cores  for 
precipitation  reconstruction  from  white  oak  at  Giant  City 
State  Park,  Illinois,  which  is  about  20  km  from  the  Dillon 
Bend  site.  Resulting  ring-with  indices  from  these  two 
studies,  available  from  the  international  Tree-Ring  Data  Bank 
(ITRDB) ,  were  plotted  against  the  Dillon  Bend  high  oak 
chronology  (Figure  12) .  The  major  wet/dry  trends  in  the 
shortleaf  pine  chronology  are  similar  to  that  of  Dillon  bend 
oak,  although  there  are  dramatic  departures  in  amplitude  for 
some  years.  A  one-year  lag  appears  in  the  final  10  years  of 
the  shortleaf  pine  chronology,  which  may  be  due  to  species - 
related  differences  or  a  missing  ring  in  the  pine  chronology, 
which  often  occurs  with  moisture  stressed  conifers  (Fritts, 
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Pigiirta  12.  Standardized  chronology  from  the  Dillon  Bend  high 
site  (several  species  of  oak)  plotted  against  TOP:  Estes 
Pine  Hills  site  (shortleaf  pine),  and  BOTTOM:  Duvick's  Giant 
City  State  Park  site. 


1976) .  Some  agreement  between  the  oak  chronologies  of  Giant 
City  and  Dillon  Bend  is  evident,  indicating  a  reflection  of  a 
common  climatic  signal. 

An  appearance  of  correlation  between  separate  study 
sites  is  encouraging,  but  most  critical  to  tree- ring  studies 
is  the  measure  of  correlation  between  growth  indices  and 
environmental  variables. 

4.B  Correlation  with  environmental  variables 

To  determine  which  single  environmental  factor  has  the 
most  influence  on  tree  growth,  it  was  necessary  to  correlate 
growth  with  all  available  environmental  measurements. 
Available  data  included  temperature,  precipitation,  and 
aggregate  measures  of  the  environment  such  as  PDSI  and  PHDI 
(see  Appendix  a) ,  and  discharge.  Initial  correlations  were 
made  with  annual  ring  indices  on  monthly  environmental  data. 
Correlations  were  made  for  the  12  months  of  the  growth  year, 
and  September  through  December  of  the  previous  year,  since 
it  is  likely  that  tree  growth  responds  to  environmental 
conditions  from  a  number  of  months,  consecutive  months  with 
relatively  high  correlation  values  were  combined  into 
seasonal  variables.  Overall,  the  highest  correlations  were 
found  between  growth  and  discharge  from  the  Mississippi 
River,  recorded  at  Thebes,  Illinois.  Further  discussion  on 
relationships  between  growth  and  environmental  variables 
follows.  Appendix  c  contains  a  summary  of  all  correlation 
results . 

4.B.1 _ Correlation  with  temperature 

Pearson  correlation  coefficients  were  computed  for  the 
floodplain,  mid- slope,  and  high-level  chronologies  to 
determine  their  relationship  with  (l)  regional  (Illinois 
Climate  Division  8) ,  (2)  Anna,  and  (3)  Carbondale  mean 
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monthly  maximum  temperature .  Mean  maximum  temperature  was 
used  since  extremely  warm  weather  is  a  limiting  factor  to 
ring  growth  during  the  growing  season.  As  expected,  computed 
relationships  were  significant  at  the  0.05  level  only  during 
the  growing  season  months.  With  all  cases,  the  high-level 
upslope  trees  were  most  sensitive  to  maximum  temperature. 

This  is  probably  due  to  the  fact  that  extremely  warm 
temperatures  during  the  growing  season  increase  evaporation 
and  moisture  stress  to  typically  moisture  sensitive  trees 
(because  of  slope  drainage) .  Therefore,  ring  growth  should 
be,  and  was  negatively  correlated  with  maximum  temperature. 

The  strongest  negative  temperature  correlations  in  this 
study  were  with  data  from  Anna  during  May  through  June  of  the 
growth  year,  with  coefficients  for  mid- slope  and  high  trees 
of  -0.31  and  -0.51,  respectively.  Results  were  not 
significant  for  the  floodplain  chronology. 

High  temperatures  in  the  latter  part  of  the  growing 
season  were  less  effective  since  the  majority  of  the  annual 
ring  forms  early  in  the  season,  while  growth  in  July  slows 
considerably  (latewood)  and  usually  ends  sometime  in  August. 
Also,  temperatures  in  the  early  spring  months  had  a  less 
significant  impact  on  growth.  These  seasonal  characteristics 
are  confirmed  by  a  dramatic  decrease  in  correlation  value  in 
the  months  prior  to  May  and  after  June  (see  Appendix  tables 
C.4  through  C.6) . 

Data  from  Anna  provided  slightly  higher  correlations 
than  Carbondale,  despite  the  fact  that  Carbondale  is  slightly 
closer  than  Anna  (20  vs.  25  km).  This  is  most  likely  due  to 
similarities  in  topography  between  Anna  and  the  Dillon  Bend 
area.  Regional  temperature  data  did  not  correlate  as  well  as 
Anna  data.  Since  temperature  tends  to  be  uniformly 
distributed  with  latitude,  the  accuracy  of  estimates  for  a 
specific  point  is  very  dependent  upon  its  location  within  the 
region.  Only  in  rare  cases  does  a  local  temperature  conform 
to  a  regional  average. 
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4.B.2  Correlation  with  precipitation 

The  spatial  distribution  of  precipitation  is  conplex, 
and  horizontally  discontinuous  (Critchfield,  1983) . 

Therefore,  it  would  be  expected  that  the  study  site,  in  the 
long  run,  would  realize  same  precipitation  characteristics 
with  Carbondale,  Anna,  and  the  region  (Illinois  Climate 
Division  8) .  This  is  indeed  the  case,  as  there  are 
significant  correlations  with  all  three  variables.  The  high 
level  chronology  had  a  0.51  correlation  coefficient  with  May- 
August  precipitation  at  Carbondale.  The  mid- slope 
chronology  best  correlated  with  Anna  May- June  precipitation 
at  0.45.  The  low  floodplain  trees  were  much  less  sensitive, 
with  no  single  dominant  season.  Floo<%>lain  correlations 
ranged  from  0.30  to  0.35  for  relationships  to  growing  season 
precipitation  and  also  a  lagged  response  to  rains  from  the 
previous  season.  Mid- slope  and  high  level  trees  also 
exhibited  correlation  to  last  growing  season's  precipitation. 
The  previous  August  consistently  showed  the  strongest  lagged 
influence,  with  correlation  coefficients  as  high  as  0.45  for 
high  level  trees. 

_ Correlation  with  EDS!  and  JBHDI 

Overall,  monthly  indices  of  PDSI  and  PHDI  for  Illinois 
Climate  Division  8  had  a  weaker  relationship  to  tree  growth 
than  precipitation,  although  there  were  significant 
correlations  at  the  0.05  level  of  confidence  for  months 
during  and  prior  to  the  growing  season.  Weaker  correlations 
are  not  surprising  since  PDSI  and  PHDI  have  an  appreciable 
lag  time  to  soil  moisture  inputs  (see  Appendix  A)  .  The 
floodplain  and  mid- slope  level  trees  had  weak  correlations 
for  all  months,  as  only  the  most  moisture  sensitive  trees 
(high  level)  had  correlation  coefficients  exceeding  0.40. 
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The  highest  correlations  with  PDSI  occurred  in  May  (0.46)  and 
June  (0.47)  of  the  growing  season  with  the  high  level  trees. 

There  were  some  relatively  strong  correlations  for  high 
level  tree  growth  and  monthly  PHDI,  but  interpretation  of  the 
results  was  difficult.  The  highest  coefficients  were  found 
with  the  current  year's  September  (0.46),  October  (0.49),  and 
November  (0.46) .  Growth  may  continue  through  September,  but 
only  rarely  continues  into  early  October  (Kramer  and 
Kozlowski,  1960) .  These  correlations  may  be  spurious  unless 
growing  season  droughts  for  this  area  naturally  continue 
through  the  fall  months. 

4 .  B_.  4.  Correlation  with  discharge 

Tree-ring  growth  indices  from  the  floodplain,  mid-slope, 
and  high  level  chronologies  were  compared  to  mean  monthly 
discharge  data  from  stream  gages  at  Murphysboro  and  Thebes, 
Illinois.  Murphysboro  is  located  10  km  upstream  from  Dillon 
Bend  on  the  Big  Muddy  River,  and  Thebes  is  51  lan  downstream 
on  the  Mississippi  River.  Correlations  between  growth  and 
Thebes  discharge  proved  to  be  much  stronger  than  those  with 
Murphysboro,  and  best  overall  from  the  tested  environmental 
variables.  Despite  the  fact  that  the  study  site  is  located 
on  the  Big  Muddy  River,  it  is  not  surprising  that  trees  there 
are  affected  by  flow  from  the  Mississippi.  As  discussed  in 
section  2.E,  the  lower  section  of  the  Big  Muddy  is  known  to 
backflood  with  high  flows  of  the  Mississippi  River. 

The  highest  Thebes  correlation  (0.57)  was  with  the  high 
level  trees  for  the  period  of  January  through  July.  Also 
very  high  was  the  floodplain  correlation  (0.53)  for  the 
period  of  previous  October  through  current  May,  making  Thebes 
discharge  the  only  variable  that  had  a  strong  relationship 
with  floodplain  trees.  The  correlation  coefficient  for  mid¬ 
slope  trees  was  not  as  impressive  (0.40)  as  results  for 
floodplain  and  high  level  trees.  As  this  was  the  case  with 
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correlations  with  other  variables,  I  suspect  that  the  small 
sample  size  of  mid- slope  trees  (9  vs.  12  floodplain  and  17 
high-level)  prevented  determination  of  a  clear  climatic 
signal. 

The  fact  that  high  level  trees  have  such  a  strong 
correlation  with  discharge  can  be  attributed  to  either  or 
both  of  two  reasons.  The  first  explanation  is  that  these 
trees  are  responding  only  to  precipitation,  which  is  the 
primary  water-balance  input  leading  to  runoff.  Second, 
despite  their  above- floodplain  location,  the  root  systems  of 
these  trees  are  sensitive  to  the  rising  and  falling  of 
floodwaters  at  Dillon  Bend.  Logic  leads  to  the  conclusion 
that  the  latter  explanation  applies.  Response  to 
precipitation  alone  is  unlikely,  because  the  discharge  at 
Thebes  results  from  runoff  occurring  in  the  entire  Upper 
Mississippi  River  watershed.  The  precipitation  (and  runoff) 
of  the  Big  Muddy  River  watershed  is  not  representative  of 
precipitation  in  the  vast  area  of  the  Upper  Mississippi 
watershed.  In  addition,  a  significant  portion  of  Mississippi 
River  runoff  is  attributed  to  spring  melting  of  the  northern 
and  alpine  snowpack.  Physical  verification  of  the  depth  of 
Dillon  Bend  root  systems  is  impossible,  but  there  are 
indications  that  they  extend  deep  enough  to  be  affected  by 
the  rise  and  fall  of  the  river.  Oak  trees  have  the 
characteristic  of  developing  a  very  long  taproot,  especially 
on  moisture-poor  sites  (Estes,  1970;  Holch,  1931;  Kramer  and 
Kozlowski,  1960) .  The  roots  of  bur  oak  have  been  measured  at 
depths  of  over  2  m  after  just  two  growing  seasons  (Holch, 
1931) ,  and  roots  of  more  mature  trees  are  known  to  extend  to 
depths  of  greater  than  10  m  (Kramer  and  Kozlowski,  1960). 

The  high  level  oak  trees  are  located  between  6  and  18  m  above 
the  lowest  floodplain  level,  with  most  trees  in  the  6-12  m 
range.  During  a  January  1993  trip  to  the  study  site, 
watermarks  were  observed  on  the  trees  and  measured  to  be  l.l 
m  above  the  low  floodplain,  most  likely  within  reach  of  at 
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least  some  of  the  high  level  tree's  roots.  That  January  1993 
flood  was  associated  with  a  stage  of  9.7  m  recorded  at 
Thebes,  not  a  major  flood  relative  to  other  known  high  flows. 
So  it  is  very  probable  that  major  backwater  floods  attained 
stages  at  Dillon  Bend  that  provided  moisture  to  the  roots  of 
mid- slope  and  high  level  trees. 

Now  that  the  best  correlated  environmental  variable 
(discharge  at  Thebes)  has  been  identified,  it  is  possible  to 
reconstruct  past  discharge  using  regression  analysis. 

4.C  Regression  analysis  results 

4..-C.1  Model  Development 

Model  development  began  with  data  retrieval  on 
dependent  and  independent  variables.  Dependent  variables 
were  the  floodplain  and  high  level  tree- ring  chronologies. 
Mid- slope  trees  were  not  included  in  regression  analysis 
because  of  weaker  correlation  with  discharge  and  other 
environmental  variables.  Independent  variables  were  Thebes 
October  (prior  year)  through  May  (growth  year)  discharge  for 
the  floodplain  trees,  and  Thebes  January- July  (growth  year) 
discharge  for  the  high  level  trees. 

The  discharge  data  set  (1934-1991)  was  divided  into  2 
equal  parts:  calibration  and  verification.  Regression 
analysis  was  conducted  using  the  Statistical  Analysis  System 
(SAS)  computer  package.  The  SAS  program  performs  an  Analysis 
of  Variance  (ANOVA) ,  conputes  adjusted  coefficient  of 
determination  (r2) ,  and  gives  regression  parameter  estimates, 
in  developing  each  regression,  output  statistics  were 
examined  for  individual  outlier  years  that  might 
disproportionately  influence  the  form  of  the  fitted 
regression  line. 


58 


Outlier  identification: 


Outliers  were  identified  if  any  of  the  following 
criteria  were  met  (recommended  by  Laboratory  of  Tree- Ring 
Research  personnel):  Cook’s  D  exceeds  0.10,  observed 
discharge  differs  more  than  2.0  standard  deviations  from 
predicted,  or  Student ized  Residual  greater  than  or  equal  to 
the  absolute  value  of  2.0. 

The  Cook's  D  test  measures  the  changes  to  the  parameter 
estimates  that  would  result  from  deleting  each  observation 
one  at  a  time,  and  Studentized  Residual  is  the  ratio  of  the 
residual  (actual  minus  predicted)  over  its  standard  error 
(Joyner,  1985) .  According  to  these  criteria,  standardized 
growth  indices  for  1943,  1968,  1975,  and  1977  were  identified 
as  outliers  in  the  floodplain  regression.  Outlier  years  in 
the  high-level  regression  were  1955,  1962,  1977  and  1983 
Investigation  of  data  for  Thebes  showed  that  there  was 
flooding  in  scane  of  the  outlier  years:  1943,  1962,  197  5,  and 
1983.  The  ring  width  indices  for  these  outlier  flood  years 
were  relatively  narrow  canpared  to  other  rings  in  the 
chronology.  Flooding  may  have  inundated  the  floodplain  to 
the  point  that  tree  growth  was  inhibited  for  those  years. 

The  other  outlier  years  (1955,  1968  and  1977)  had  relatively 
low  flow  all  through  the  growing  season,  but  yet  fairly  wide 
growth  rings.  Carbondale  climate  data  for  those  years  shows 
that  precipitation  was  near  long-term  average  during  the 
growing  season,  so  trees  probably  did  not  experience  moisture 
stress,  and  realized  healthy  ring  growth. 

Parameter  estimates: 

Each  of  these  outlier  years  were  eliminated  from  the 
data,  and  regression  equations  were  finalized.  Regression 
equations  for  the  two  halves  of  the  data  set,  (1934-1963)  and 
(1964-1991),  were  evaluated  by  their  values  of  explained 
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variation,  or  adjusted  r2.  The  early  period  (1934-1963) 
exhibited  slightly  higher  values  of  adjusted  r2  for  both  the 
floodplain  (0.41)  and  high  level  (0.45)  chronologies,  so  it 
was  used  to  calibrate  the  model.  The  later  period  data 
(1964-1991)  was  set  aside  for  use  in  verification  of  the 
model.  The  final  regression  equations  to  estimate  October- 
May  and  January- July  Thebes  discharge  amounts  were 

Q (October-May)  =  -24395  +  67352X  (1) 

Q( January -July)  =  -17404  +  59664X  (2) 

where  X  is  the  ring  width  index  for  a  particular  year,  and  Q 
is  total  discharge  (m3/s)  for  the  seasonal  period.  Values  of 
the  F-ratio  for  (l)  and  (2)  are  18.9  and  21.3,  respectively. 
Their  acceptance  at  the  99.9%  level  of  confidence  assures 
that  the  slope  parameters  are  signif icantly  different  from 
0.0.  It  was  next  necessary  to  test  the  accuracy  of  (l)  and 
(2)  as  predictors  of  seasonal  discharge. 

4.C.2  Verification  of  regression  model 

To  test  the  accuracy  of  the  regression  equations,  each 
equation  was  used  to  predict  seasonal  Thebes  discharge  for 
the  1964-1991  period.  Using  the  computer  program  VERIFY 
(Holmes,  et  al.,  1986),  the  predicted  values  of  discharge 
were  then  compared  with  actual  data.  The  floodplain  and 
high-level  regression  equations  passed  all  verification 
criteria,  as  listed  in  Table  ll. 

Using  these  equations,  it  was  possible  to  reconstruct 
seasonal  discharge  at  Thebes  through  the  length  of  the 
floodplain  or  upland  tree- ring  chronologies.  After  extending 
the  hydrologic  record,  patterns  and  trends  of  discharge 
variability  at  at  Thebes  were  examined.  To  do  this, 
reconstructed  and  actual  Thebes  discharge  was  plotted  using 
both  the  high  and  low  chronology  regression  equations.  The 
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Verification  of  Floodplain  Regression 


VALUE 

OBTAINED 


Correlation  0.48 
Reduction  of  error  0.24 
T- value  1.80 
Sign-products  6 


VALUE  OF 
95%  CONFIDENCE 

£  0.32  passed 

£0.10  passed 

£  1.71  passed 

£  9  passed 


Verification  of  High  Level  Regression 


VALUE 

OBTAINED 

Correlation  0 . 56 

Reduction  of  error  0.39 
T- value  2.16 

Sign-products  6 


VALUE  OF 
95%  CONFIDENCE 

£0.32  passed 

£  0.10  passed 

£  1.71  passed 

1  9  passed 


Tahi ft  i i -  Results  of  verification  of  the  floodplain  and  high 
level  regression  equations,  based  on  comparison  of  confuted 
test  statistics  and  verification  criteria. 
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Discharge  at  Thebes,  Illinois 


Reconstructed - Actual 


pigure  13 .  Actual  and  reconstructed  January  through  July 
discharge  for  Thebes,  Illinois,  from  regression  based  on 
high-level  trees. 
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reconstruction  from  high  level  trees  provided  a  more 
extensive  record,  dating  to  1869,  and  was  the  focus  of 
analysis.  The  analysis  began  with  comparing  discharge  plots 
from  the  gaged  period  (1934-1991) (Figure  13) . 

Trends  in  the  gaged  portion  of  the  plot  are  in  general 
agreement,  but  the  amplitude  of  reconstructed  discharge  is 
often  less  than  actual  extremes  of  flow.  Characteristics  of 
regression  analysis  provide  a  plausible  explanation.  Values 
of  adjusted  R2  for  the  floodplain  and  high  level  regressions 
were  41%  and  45%,  respectively.  As  values  of  R2  become 
small,  the  variance  explained  by  the  regression  decreases. 
This  results  in  a  reduction  in  amplitude  in  reconstructed 
data  versus  actual  data. 

The  reconstructed  plot  exhibits  a  lagged  response  to  the 
regional  droughts  of  1980-81  and  1988,  but  responds  quickly 
to  the  increase  in  flows  during  the  wet  years  of  1982-86 
(Changnon,  1991) .  It  must  be  recognized  that  the  extreme 
fluctuations  in  the  earliest  reconstructed  decade  (187 Os)  may 
be  less  reliable  than  later  years,  since  only  l  tree 
contributed  to  the  chronology  for  that  time. 

Discharge  variability  in  the  1900s  was  marked  by 
decadal -scale  trends,  but  overall  the  range  was  limited 
between  20k  and  90k  m3/s.  The  years  1870-1900  contain  the 
highest  and  lowest  reconstructed  flows  for  the  entire  period. 
Inclusion  of  flow  amplitude  data  such  as  these  could  have  an 
effect  on  flood- frequency  computations,  and  possibly  call  for 
a  re-evaluation  of  flood  control  projects.  To  investigate 
this  possibility,  the  actual  and  reconstructed  January 
through  July  discharge  were  used  to  construct  flow  frequency 
curves  (Figure  14) .  in  general,  the  reconstructed 
probabilities  slightly  underestimate  the  magnitudes  of  the 
actual  record.  An  exception  occurs  in  the  region  of  low 
probability  events,  inclusion  of  higher  reconstructed  flows 
significantly  increases  the  magnitude  of  rare  events  over 
what  a  linear  extension  of  the  actual  flow  curve  would  have 
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ow  frequency  plots  constructed  of  actual  (marked  by  x)  and  reconstructed 
January  through  July  discharge  for  Thebes,  Illinois. 


estimated.  From  this  approach,  the  actual  probability  curve 
may  be  underestimating  the  magnitude  of  rare  events  such  as 
the  100  year  or  greater  floods. 

Another  approach  locates  the  high  reconstructed  flows  on 
an  extrapolation  of  the  reconstructed  curve.  This  action 
would  place  them  as  400  to  600  year  events.  Although  neither 
approach  can  be  proven,  the  second  method  is  preferred  since 
probability  curves  are  assumed  to  maintain  a  smooth  form 
without  abrupt  changes  (Kite,  1977).  With  this 
consideration,  implications  of  including  reconstructed  flows 
in  flood  control  project  design  would  be  less  severe. 

A  third  approach  concludes  the  high  reconstructed  flows 
of  the  late  1800s  to  have  been  part  of  an  abnormally  high 
flow  period  (i.e.  1850-1900) .  Flow  frequency  based  on  data 
from  the  abnormally  high  period  may  be  very  different  from 
frequency  from  the  gaged  period  (1934 -present) .  The  best 
conclusion  realizes  that  there  may  have  been  increased 
variability  of  prehistoric  flows  over  what  has  been  observed 
in  the  gaged  period.  The  abnormally  high  flow  period  must  be 
confirmed  through  further  tree  ring  or  other  proxy  analysis. 

Confidence  in  the  discharge  reconstruction  model  earn 
also  be  gained  by  confirming  pre-historic  floods  with  other 
data  sources.  Flood- ring  analysis  at  the  Dillon  Bend  site 
provides  evidence  for  a  number  of  significant  flood  events. 

4 . D  Flood  ring  analysis 

All  of  the  floodplain  trees  (oak  and  hackberry)  and  mid¬ 
slope  trees  (oak)  were  carefully  inspected  under  a  microscope 
for  flood  rings  (see  section  l.D.l).  The  entire  length  of 
all  cores  was  examined,  but  rings  near  the  pith  and  bark  were 
too  narrow  to  conclusively  locate  flood  rings.  As  a  result, 
no  dated  flood  rings  precede  the  start  of  stage  gage 
operation  (1896) .  Therefore,  the  main  purpose  of  the 
analysis  was  to  determine  the  effectiveness  of  the  flood-ring 
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reconstruction  technique,  and  its  possible  application  to 
specific  locations  on  floodplains,  such  as  the  Dillon  Bend 
section  of  the  Big  Muddy  River. 

From  the  examined  cores,  growth  years  of  discovered 
flood  rings  were  used  to  prepare  a  flood  frequency  chart 
(Table  12,  Figure  15).  There  were  83  rings  with  signs  of 
flooding,  all  from  floodplain  trees.  Cores  reflected  the 
floods  of  1947  and  1951  most  often,  with  rings  for  that  year 
exhibiting  a  distinct  band  of  enlarged  vessels  near  the  end 
of  latewood  growth.  The  intra-ring  position  of  abnormal 
cells  corresponded  roughly  to  timing  of  the  flood  within  the 
growing  season,  considering  a  growing  season  that  begins  in 
March  and  ends  in  August.  For  example,  an  flood  ring  located 
in  the  earlywood  zone  may  be  related  to  a  March  or  April 
flood  (Table  12) .  Such  interpretations  are  loose  and  must  be 
based  on  growing  season  length  for  a  particular  site. 

The  83  flood  rings  corresponded  to  25  individual  growth 
years,  each  of  which  corresponded  to  an  annual  peak  stage  of 
at  least  10.3  m  documented  at  Thebes  (Table  12);  slightly 
greater  than  the  9.7  m  "rule  of  thumb"  for  overbank  flooding 
on  the  lower  Big  Muddy  River  (Ballance,  pers.  comm.,  1992). 
Gaged  records  reveal  that  the  stage  at  Thebes  equalled  or 
exceeded  9.7  m  in  45  years  over  the  stage  recorder's  period 
of  operation  (1896-1991),  with  41  of  the  floods  occurring 
during  the  growing  season  (approximately  March  through 
September) .  There  are  a  number  of  possible  reasons  why  only 
61%  of  the  growing  season  flood  years  were  detected. 

First  of  all,  response  of  trees  are  likely  dependent 
upon  flood  magnitude.  It  is  unlikely  that  flood  magnitude 
was  great  enough  inundate  tree  crowns,  unless  the  flooding 
occurred  early  in  the  tree's  life.  So,  defoliation  flood 
rings  did  not  exist.  The  source  of  flood  evidence  is  left  to 
white  flood  rings  caused  by  longer- duration,  but  non¬ 
inundating)  events.  The  key  factor  is  a  flood  duration  long 
enough  to  cause  rapid  shoot  growth  (Yanosky,  1983) . 
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Figure  15.  Time  scale  plot  of  f 
examining  floodplain  oak  and  hac 
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Tahi <3  i ?.  Information  on  flood  rings  from  floodplain  trees 
at  Dillon  Bend,  included  are:  year,  intra-ring  position, 
number  of  verifying  cores,  and  stage  and  date  of  related 
flooding  at  Thebes.  Codes  for  intra-ring  position  are:  EW 
(earlywood) ,  ELW  (early- latewood) ,  MLW  (mid- latewood) ,  LLW 
(late-latewood) . 
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The  final  phase  of  flood  ring  analysis  was  stage 
investigation.  Since  the  floodplain  at  Dillon  Bend  has  about 
l  m  of  relief  relative  to  the  river  surface,  a  small 
magnitude  event  would  most  likely  produce  flood  rings  only  in 
the  lowest  trees.  A  higher  magnitude  event  may  produce  flood 
rings  in  both  elevated  and  low  trees.  The  core  which  held 
the  most  flood  rings  was  from  the  lowest -lying  tree 
(hackberry) .  The  lowest  trees  had  the  benefit  of  maximized 
duration,  being  affected  through  the  rising  and  falling  of 
the  floodcrest  while  higher  trees  may  be  flooded  only 
briefly.  Comparison  between  tree  elevation  and  magnitude 
data  did  not  reveal  any  clear  stage  influence,  probably 
because  the  duration  of  floodwaters  at  some  threshold  level 
is  most  inportant  for  flood  ring  formation.  It  was 
interesting  to  note,  however,  that  there  was  a  distinct  level 
above  which  trees  no  longer  contained  flood  rings.  This 
began  at  a  point  2.7  m  above  the  lowest  level  of  the 
floodplain.  It  can  be  concluded  that  floods  of  damaging 
duration  have  not  risen  above  that  stage  over  the  period  of 
the  tree- ring  chronology  at  this  site. 

Flood  ring  information  provides  a  means  of  relating 
gaged  stage  data  to  resulting  water  level  at  other  locations, 
or  reconstructing  flood  history  for  a  river  lacking 
hydrologic  data.  There  are  obvious  benefits  from  this  work 
for  agriculture  and  construction,  for  those  industries  are 
highly  dependent  on  flood  stages  in  fields  or  developed 
areas,  not  on  the  actual  discharge  at  some  distant  gaging 
station.  The  flood  ring  analysis  technique  is  especially 
useful  on  rivers,  such  as  the  Big  Muddy,  that  are  influenced 
by  flow  from  another  drainage  basin. 

4.B  Results  of  elevational  comparisons 

The  initial  step  in  an  elevational  analysis  of  ring 
widths  was  computation  of  correlation  coefficients  between 
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floodplain/high-level  index  ratios  (dependent  variable)  and 
hydrologic  data  from  Thebes,  Illinois.  Peak  annual  stage  was 
chosen  as  the  independent  variable,  since  the  maximum  height 
of  floodwaters  determines  degree  of  tree  inundation  more  than 
a  seasonal  aggregate  of  discharge.  Results  in  the  linear 
correlation  analysis  were  not  significant  at  the  95%  level  of 
confidence,  and  it  was  evident  that  an  alternate  approach  to 
the  problem  was  necessary. 

when  considering  the  floodplain/high- level  ratio 
concept,  it  is  important  to  recognize  the  possibility  of 
thresholds.  That  is  to  say,  there  may  be  a  threshold  of  peak 
stage  above  which  floodwaters  will  cause  inhibited  ring 
growth  in  floodplain  trees.  Also  there  may  be  a  threshold  of 
the  ratio  itself;  a  significant  departure  from  1.0  which 
corresponds  to  abnormal  comparative  growth  from  flooding. 

The  threshold  concept  can  be  investigated  through  the  use  of 
categorical  variables  and  Analysis  of  Variance  (ANOVA) ,  which 
will  conclude  whether  or  not  there  are  explained  differences 
between  the  categories.  The  ANOVA  approach,  and  its 
underlying  assumptions  will  be  explained  in  detail  in  the 
remainder  of  this  section. 

_ Assumptions  of  .ANQ va 

The  One-Way  ANOVA  test  is  concerned  with  differences 
between  classes,  and  uses  the  class  means,  sums  of  squares, 
and  estimates  of  variance  to  summarize  their  characteristics 
(Clark  and  Hosking,  1986) .  Total  variance  in  the  data  is 
divided  into  explained  (between  class)  and  unexplained  parts 
(within  class) .  A  ratio  (F)  of  between/within  class  variance 
is  determined.  If  between- class  variance  is  significantly 
greater  than  within- class  variance,  the  null  hypothesis  of 
"no  difference  between  categories"  will  be  rejected.  A  level 
of  confidence  of  95%  was  chosen  as  a  reasonable  interval 
estimate  for  parameters. 
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The  two  variables  used  in  the  ANOVA  were  ratio,  and 
stage.  The  main  principle  behind  using  the  floodp lain/high- 
level  ratio  is  that  values  of  ratio  near  1.0  indicate  similar 
growing  conditions  between  the  two  sites.  Positive 
departures  from  1.0  indicate  enhanced  floodplain  growth  from 
flood-delivered  moisture,  while  negative  departures  show 
inhibited  floodplain  growth  from  inundation. 

The  degree  of  departure  required  for  each  of  these 
effects  was  estimated  by  adding  and  subtracting  the  value  of 
0.10,  which  is  10%  of  the  equal- site  growth  ratio,  to  the 
equal- site  growth  ratio  (1.0) .  The  resulting  threshold  for 
flood  enhanced  growth  was:  ratio  £  1.10,  and  for  flood- 
inhibited  growth:  ratio  £  0.90.  When  ratio  meets  either  of 
these  criteria,  flooding  is  determined  to  have  occurred  in 
that  year.  Threshold  for  the  peak  stage  variable  was  set  at 
9.7  m,  since  the  river  is  known  to  overflow  the  banks  when 
the  Thebes  stage  reaches  that  level  (Ballance,  pers.  comm., 
1992) .  All  years  with  stages  greater  than  9.7  m  were  labeled 
flood  years. 

Each  year  of  the  chronology  was  listed  with  its 
corresponding  value  of  ratio  and  stage.  The  variables  were 
then  assigned  categorical  variables  of  either  l  (flooding)  or 
0  (no  flooding)  according  to  the  criteria  previously 
described  (see  Table  13) .  When  categorical  values  agree  for 
a  year  (both  0  or  l) ,  the  ratio  measure  accurately  recorded 
the  presence  or  absence  of  a  flood.  With  the  thresholds 
determined  for  Dillon  Bend,  the  variables  agreed  66%  of  the 
time,  and  70%  Thebes  floods  at  or  above  9.7  m  were  detected. 

Most  of  the  errors  (13  of  19)  occurred  when  an  actual 
flood  was  not  reflected  by  the  ratio  variable.  The  probable 
cause  is  that  those  missed  floods  were  not  of  sufficient 
duration  to  alter  ring  growth.  There  were  4  cases  of  small 
ratio  values  with  no  flooding.  It  is  possible  that  dry 
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1934  0.88  1  6.2  0  No  1970 
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Table  13 .  values  of  categorical  variables  ratio  and  stage,  and  prediction  agreement 
statistics  for  each  year  of  the  combined  floodplain/high-level  chronology. 


conditions  did  not  affect  trees  on  the  slope  as  much  as  those 
on  the  floodplain,  since  slope  trees  are  more  adapted  to 
stressful  environmental  conditions.  In  only  2  cases  when  the 
ratio  was  £  1.10  did  flooding  not  occur. 

4..E.3 _ Results  of  ANOVA 

ANOVA  was  run  on  the  categorical  variables  ratio  and 
stage,  to  determine  if  there  is  a  difference  between  the  non¬ 
flooding  ratio  (0.90  <  ratio  <  1.10)  and  the  flooding  ratio 
(0.90  >  ratio  £  1.10).  The  null  hypothesis.  Ho,  is  that 
there  is  no  difference  between  the  populations  for  the  two 
ratio  classes.  ANOVA  results  are  shown  in  Table  14.  The 
rejection  region  at  the  95%  level  of  confidence  is  F  £  4.0. 
Since  the  calculated  F  is  5.76,  the  null  hypothesis  is 
rejected  and  it  is  concluded  that  there  is  a  significant 
difference  between  the  two  ratio  categories. 

The  preliminary  results  of  the  comparative  technique  are 
promising,  especially  since  it  detected  Dillon  Bend  floods 
slightly  better  than  the  flood  ring  technique  (66%  vs.  61%) . 
Both  approaches  have  potential  for  improved  success,  and 
should  be  tested  at  other  locations. 

Analysis  of  Variance 

Source  of  variation _ SS _ df  MS _ F  Prob  >  F 

Between  groups  1.33  1  1.34  5.76  0.02 

Within  groups _ i3.oi  56  0.23 _ 

Total  14.34  57 

Table  14 .  Results  of  Analysis  of  Variance  between 
categorical  measures  of  the  floodplain/high-level  ratio 
variable,  and  a  categorical  measure  of  peak  annual  stage. 
(SS=sum  of  squares;  MS=mean  sum  of  squares;,  df=degrees  of 
freedom;  F=MSfcetween/MSwithin) • 
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CHAPTER  5 

SUMMARY  AMD  CONCLUSIONS 


It  is  possible  to  reconstruct  flood  events  on  the  Big 
Muddy  River  using  any  of  3  methods:  1)  regression  of  ring 
width  with  gaged  discharge;  2)  flood  ring  analysis;  and  3) 
floodplain/upland  tree  comparison.  The  regression  approach 
models  Thebes  seasonal  (i.e.  January  through  July)  discharge 
as  a  function  of  standardized  ring-width  indices  of  upslope 
trees.  Flood  ring  analysis  associates  abnormalities  in  the 
wood  tissue  of  a  growth  ring  to  physiological  impacts  of 
flooding.  Comparison  of  floodplain  and  upland  trees  relies 
on  the  premise  that  floodplain  trees  have  inhibited  ring 
growth  when  a  flood  event  results  in  prolonged  soil 
saturation,  while  upland  trees  have  enhanced  or  normal 
growth.  Each  of  these  approaches  in  dendrohydrology  contains 
unique  information  on  past  flooding  and  holds  special 
potential  application  for  Illinois  and  other  midwestem 
states,  since  much  of  the  region's  development  is  located 
along  waterways  within  the  Mississippi  River  basin. 

Regression  of  ring  width  on  discharge  gives  an 
indication  of  the  overall  hydrology  (i.e.  low  flow,  high 
flow)  of  the  years  which  comprise  the  tree- ring  chronology. 

A  time  line  of  aggregate  discharge  is  constructed,  from  which 
patterns,  trends,  and  extremes  can  be  extracted.  Flood  ring 
analysis  and  floodplain/upland  tree  conparisons  are  mainly 
used  to  determine  a  prehistoric  flood  chronology  for  the 
study  site.  Rings  which  are  determined  to  reflect  flooding 
are  assigned  calendar  years  which  can  be  used  to  construct  a 
flood- frequency  plot  for  the  site.  Flood  ring  analysis  can 
also  provide  some  flood  magnitude  information,  by  noting  the 
elevation  at  which  trees  no  longer  exhibit  flood  rings.  When 
the  3  described  approaches  to  flood  reconstruction  are 
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applied  simultaneously,  many  characteristics  of  flooding  for 
a  particular  site  are  revealed. 

The  study  site,  located  along  the  Dillon  Bend  section  of 
the  Big  Muddy  River,  was  cored  extensively  for  flood 
reconstruction.  A  total  of  62  cores  were  collected  from  23 
oak  trees  and  8  hackberry  trees.  The  study  area  was  divided 
into  3  separate  collection  sites:  floodplain,  mid- slope  and 
high-level  trees.  Standardized  ring-width  indices  from  the 
floodplain  and  high-level  trees  were  regressed  on  seasonal 
aggregates  of  Thebes  discharge,  in  order  to  compute  a 
discharge  reconstruction  model.  The  resulting  equations, 
when  used  with  the  ring -width  data,  extend  seasonal  Thebes 
discharge  information  (January  through  July)  from  1934  to 
1870.  The  extended  discharge  data  shows  that  the  years 
1870  to  1900  contained  the  highest  and  lowest  flow  seasons 
for  the  entire  period  (1870  to  1991) .  Inclusion  of 
reconstructed  discharge  into  flood  frequency  calculations 
reveals  that  studies  based  on  gaged  data  may  '  a 
underestimating  the  magnitude  of  rare  flood  event.-,  at  Thebes. 

Inspection  of  floodplain  cores  for  flood  rings  lead  to 
the  discovery  of  83  such  abnormal  rings.  Comparison  of 
tabulated  flood  ring  results  to  Thebes  stage  data,  showed 
that  the  technique  detected  61%  of  floods  above  9.7  m.  There 
was  a  distinct  level  on  the  landscape  above  which  trees  no 
longer  contained  flood  rings,  giving  a  rough  estimate  of  a 
local  peak  stage  for  the  length  of  the  chronology  (87  years) . 

The  elevational  comparison  technique  was  more  successful 
in  detecting  flood  events  than  the  flood  ring  technique  (70% 
vs.  61%) .  A  ratio  of  floodplain/high-level  standardized 
ring -width  indices  was  calculated  and  compared  to  peak  annual 
stage  data  for  Thebes.  The  ratio  and  stage  data  were 
converted  into  categorical  variables  for  analysis  of  variance 
(ANOVA) .  It  was  determined  that  there  is  a  difference 
between  flooding  and  non- flooding  ratio  variables,  which 
detected  70%  of  the  Thebes  floods  above  9.7  m. 
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Considering  that  flooding  is  one  of  the  most  damaging 
and  costly  weather- related  disasters  in  the  United  States, 
the  importance  of  successful  floodplain  management  cannot  be 
disputed.  Since  funding  for  management  is  limited,  it  is 
necessary  to  design  measured  and  practical  flood  control 
systems.  Current  methodology  uses  stream  gage  data  to 
estimate  the  magnitude  and  frequency  of  floods  for  a 
particular  river.  When  this  method  is  based  on  a  short 
period  of  record,  which  is  the  case  for  much  of  Illinois  and 
the  United  States,  evaluation  of  the  magnitude  of  rare  flood 
events  becomes  difficult.  This  problem  can  be  alleviated 
with  the  inclusion  of  proxy  flood  data  from  other  sources. 
Tree- ring  analysis  has  proved  to  be  a  valuable  source  of 
paleoflood  data.  Inclusion  of  reconstructed  discharge  in 
flood  frequency  estimates  could  prompt  a  reassessment  of 
existing  flood  control  projects,  and  influence  the  design  of 
future  protection.  Dendrohydrology  also  has  the  flexibility 
of  giving  past -flood  information  at  a  particular  site,  which 
may  be  far  removed  from  a  gaging  station.  The  economic 
advantage  to  developers,  farmers,  and  current  residents  are 
obvious,  in  addition  to  the  potential  savings  of  lives  and 
property. 

With  this  in  mind,  dendrohydrology  should  continue  to  be 
tested  on  other  rivers,  with  other  tree  species,  and  new 
techniques  should  be  investigated. 
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Appendix  A:  PDSI  and  PHD! 


Palmer  Drought  Severity  Index  (PDSI) : 

PDSI  is  a  method  whereby  the  agricultural  effects  of  a 
meteorological  drought  cam  be  calculated  from  temperature  and 
precipitation,  in  its  calculation,  drought  severity  is 
viewed  as  a  function  of  moisture  demand  and  moisture  supply 
(see  Fieldhouse  and  Palmer,  1965) .  Therefore,  it  determines 
whether  or  not  the  water  needs  of  crops  (or  trees)  at  a 
particular  location  are  being  met.  Since  it  takes  time  for 
serious  drought  to  develop,  PDSI  takes  into  account  previous 
and  current  weather  conditions.  The  developed  procedure 
treats  drought  severity  as  a  function  of  accumulated 
differences  between  actual  and  required  precipitation,  where 
the  requirement  depends  on  rainfall  carryovers,  evapo- 
transpiration,  moisture  recharge,  and  runoff  (Fieldhouse  and 
Palmer,  1965)  .  The  computed  monthly  drought  severity  index 
(X)  range  from  about  -5.0  to  +5.0,  and  is  interpreted  by  its 
class  of  wet  and  dry  periods  (Table  A.l) . 


Monthly  value,  of  X  Clans 


£4.00 

extremely  wet 

3.00 

to  3.99 

very  wet 

2.00 

to  2.99 

moderately  wet 

1.00 

to  1.99 

slightly  wet 

0.50 

to  0.99 

incipient  wet  spell 

0.49 

to  -0.49 

near  normal 

-0.50 

to  *0.99 

incipient  drought 

-i.eo 

to  -1.99 

mild  drought 

-2.00 

to  -2.99 

moderate  drought 

-3.00 

to  *3.99 

severe  drought 

S  -4.00 

extreme  drought 

Table  A.l  PDSI  classes  for  wet  and  dry  periods. 
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Palmer  Hydrologic  Drought  Index  (PHDI) : 


The  "moisture  anomaly  index"  (PHDI,  or  Z- index)  is 
sensitive  to  short-term  changes  in  soil  moisture,  and 
responds  more  quickly  than  PDSI.  Therefore,  PHDI  is  a  better 
indicator  of  agricultural  drought  than  PDSI,  since  crop 
growth  is  highly  dependent  upon  current  moisture  levels 
during  the  growing  season  (Isard  and  Easterling,  1989) .  Tree 
growth  is  most  sensitive  to  surface  moisture  supply,  but 
because  of  extensive  root  systems,  can  also  be  affected  by  a 
lowering  of  the  water  table.  As  a  result,  tree  growth  may 
reflect  long  term  changes  in  the  water  table  as  measured  by 
PDSI. 
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Appendix  B:  Computer  programs 


Computer  programs  used  in  this  research  were  obtained 
from  the  Laboratory  of  Tree- Ring  Research,  University  of 
Arizona,  Tucson.  They  were  developed  by  Richard  Holmes  and 
Edward  Cook  for  the  international  Tree-Ring  Data  Bank 
(ITRDB) .  A  short  description  of  each  program  follows: 

COFECHA: 

The  main  function  of  COFECHA  is  identification  of  dating 
errors  in  tree- ring  measurements  (Holmes,  1983) .  The  dated 
and  measured  ring  series  are  filtered  by  fitting  a  20 -year 
cubic  spline,  and  then  dividing  the  series  values  by  the 
corresponding  spline  curve  values  to  remove  low-frequency 
variance.  The  high  frequency  residual  is  subjected  to  a  log 
transformation  to  equalize  proportionally  the  variability 
among  small  and  large  rings  (Holmes,  et  al.,  1986).  The  mean 
values  of  the  combined  series  comprise  the  master  dating 
chronology.  After  an  initial  COFECHA  run  using  all  of  the 
ring  series,  the  10  series  with  the  best  correlation  to  the 
master  are  selected  to  compose  a  more  reliable  master  dating 
chronology.  Each  COFECHA  output  identifies  measuring  or 
counting  errors  by  flagging  40  year  ring  segments  and 
individual  outliers  within  each  core  that  do  not  agree  with 
the  master.  After  using  the  output  to  account  for  possible 
problems  in  master  cores,  each  of  the  remaining  cores  can  be 
added  one  by  one  to  the  COFECHA  program  and  checked  against 
the  master  chronology.  As  measuring  or  counting  errors  are 
remedied,  each  core  is  included  in  the  final  dating 
chronology. 
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ARSTAN: 


The  program  ARSTAN  produces  annual  standardized, 
dimensionless  indices  for  all  rings  of  one  core.  The  first 
step  is  to  fit  either  negative  exponential  curves  or  stiff 
cubic  splines  to  tree-ring  series  (Holmes,  et  al.,  1986). 
Trend  curves  differ  in  form  depending  upon  site- specific 
growth  variations  and  distance  from  the  pith  of  the  tree, 
visual  inspection  of  raw  ring-widths  is  used  to  estimate 
whether  the  trend  of  growth  is  a  clear  departure  from 
negative  exponential.  ARSTAN  prompts  as  to  whether  the 
operator  believes  there  are  any  cores  with  clear  departures. 
Specified  cores  are  fitted  with  cubic  splines  to  model  growth 
trend.  An  annual  index  is  computed  by  dividing  the  measured 
ring  width  by  the  value  of  the  growth  curve  for  that 
particular  year.  To  account  for  the  lag  influence  of 
previous  years'  growth,  ARSTAN  uses  autoregressive  modeling 
to  enhance  the  common  signal  (Holmes,  et  al.,  1986).  The 
chronology  of  a  tree  core  is  then  expressed  as  a  series  of 
indices  that  are  stationary  in  both  mean  and  variance 
(Stockton,  1975) . 

Master  chronologies  are  constructed  by  combining  the 
indices  from  individual  cores,  and  compur ing  an  average 
growth  index  for  each  year  in  the  combined  chronology. 

VERIFY: 

The  purpose  of  the  VERIFY  program  is  to  test  the 
accuracy  of  the  discharge  reconstruction  model.  This  is 
accomplished  by  comparing  predicted  values  of  discharge  to 
actual  data  that  were  not  included  in  the  determination  of 
the  regression  parameters.  The  following  test  statistics  are 
computed  for  the  data  sets  at  the  95%  level  of  confidence: 
correlation  coefficient,  reduction  of  error,  T-value,  and 
sign -products.  If  any  of  the  computed  statistics  fail  to 


86 


meet  acceptance  criteria,  the  discharge  reconstruction  model 
fails. 

MEASURE: 

The  MEASURE  software  tests  the  accuracy  of  ring 
measurement  by  comparing  original  measurements  with 
remeasured  ones  for  the  same  core.  The  program  prompts  the 
user  for  the  ring-width  series  for  the  original  and 
remeasured  cores.  Descriptive  statistics  (mean,  standard 
deviation,  sensitivity,  etc.)  are  computed  for  each  data  set. 
A  difference  of  means  test  is  conducted  at  the  95%  and  99% 
level  of  confidence,  and  printed  results  show  whether  cr  not 
these  standards  are  met. 
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Appendix  C 


Complete  results  of  correlation  analysis 

Results  of  linear  correlation  analysis  between 
floodplain,  mid- slope,  and  high-level  growth  indices  and 
environmental  variables  are  presented  in  tabular  format 
(Tables  C.l  through  C.10) .  Presentation  of  all  correlation 
results  are  important  to  show  that  tree  rings  do/do  not 
respond  to  seasonal  environmental  variables  in  explainable 
degree  and  fashion.  Each  table  corresponds  to  a  different 
environmental  variable: 

Table _ Environmental  variable 

C.l  Illinois  regional  precipitation 

C. 2  Anna  precipitation 

C.3  Carbondale  precipitation 

C.4  Illinois  regional  temperature 

C . 5  Anna  temperature 

C . 6  Carbondale  temperature 

C.7  Illinois  regional  PDSI 

C. 8  Illinois  regional  PHDI 

C. 9  Murphysboro  discharge 

C. 10  Thebes  discharge 

Each  row  of  data  are  results  for  one  chronology:  LOW 
(floodplain) ,  MID  (mid- slope) ,  and  HI  (high-level).  Listed 
in  each  row  are  linear  correlation  coefficients  (r) ,  along 
with  the  probability  of  the  null  hypothesis  (true  population 
correlation  coefficient  *  0.0),  below  the  correlation 
coefficient.  A  95%  level  of  confidence  corresponds  to  a 
probability  £  0.05. 

Months  analyzed  are  October  through  December  after  the 
past  growth  season,  and  January  through  September  of  the 
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Illinois  regional  precipitation 

Pearson  Correlation  Coefficients  /  Prob  >  |t|  under  Mo:  »ho*0  /  Nustber  of  Observations 
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Tahiti  r.i  Correlation  coefficients,  and  their  level  of  significance,  for  standardized 
ring  width  and  Illinois  regional  precipitation  (climate  division  8) . 

Low  =  floodplain;  mid  =  mid- slope;  hi  =  high-level; 
prefix  L  -  last  year. 


Tablfi_£^  Correlation  coefficients,  and  their  level  of  significance,  for  standardized 
ring  width  and  Anna  precipitation. 

Low  *  floodplain;  mid  =  mid- slope;  hi  =  high-level; 
prefix  L  =  last  year. 


Carbondale  precipitation 
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Tahlft  Correlation  coefficients,  and  their  level  of  significance,  for  standardized 

ring  width  and  Carbondale  precipitation. 

Low  =  floodplain;  mid  =  mid-slope;  hi  =  high-level; 
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Tahlp  p . 5  Correlation  coefficients,  and  their  level  of  significance,  for  standardized 
ring  width  and  Anna  tenqperature . 

Low  =»  floodplain;  mid  =  mid- slope;  hi  =  high-level; 
prefix  L  =  last  year. 


Rho«0  /  Number  of  Observations 
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Correlation  coefficients,  and  their  level  of  significance,  for  standardized 
ring  width  and  Carbondale  temperature . 

Low  «  floodplain;  mid  *  mid- slope;  hi  =  high-level; 
prefix  L  -  last  year. 


Pearson  Correlation  Coefficients  /  Prob  >  |R|  under  Nos  Rho»0  /  Number  of  Observations 
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Table  c.io  (continued) 


Appendix  D 


Standardized  indices  for  chronologies 


Listed  are  ring  growth  years,  and  corresponding 
standardized  dimensionless  indices  for  the  floodplain  (LOW) , 
mid- slope  (MID),  and  high-level  (HIGH)  chronologies. 


Year  LOW  MID.  HIGH 


Year  LOW -MID  , HIGH 


Year  LOW  mid  HIGH 


1870 

1871 

1872 

1873 

1874 

1875 

1876 

1877 

1878 

1879 

1880 
1881 
1882 

1883 

1884 

1885 

1886 

1887 

1888 

1889 

1890 

1891 

1892 

1893 

1894 

1895 

1896 

1897 

1898 

1899 

1900 

1901 

1902 

1903 

1904 

1905  153 

1906  114 

1907  9S 

1908  92 

1909  57 

1910  64 


205 

1911 

57 

46 

71 

1952 

131 

103 

95 

150 

1912 

154 

73 

121 

1953 

79 

94 

87 

65 

1913 

146 

39 

79 

1954 

93 

104 

92 

60 

1914 

97 

42 

75 

1955 

109 

136 

122 

46 

1915 

140 

71 

132 

1956 

82 

108 

99 

76 

1916 

104 

70 

123 

1957 

103 

111 

99 

59 

1917 

110 

91 

115 

1958 

103 

113 

105 

90 

1918 

82 

60 

84 

1959 

82 

101 

87 

94 

1919 

72 

102 

109 

1960 

110 

108 

106 

145 

1920 

124 

92 

113 

1961 

97 

92 

86 

121 

1921 

121 

85 

102 

1962 

97 

95 

82 

113 

1922 

124 

172 

124 

1^63 

83 

84 

69 

72 

1923 

91 

146 

121 

1964 

86 

85 

74 

93 

1924 

97 

157 

146 

1965 

111 

103 

93 

117 

1925 

95 

118 

115 

1966 

88 

101 

92 

59 

1926 

91 

109 

90 

1967 

104 

105 

102 

44 

1927 

112 

175 

122 

1968 

119 

123 

109 

43 

1928 

106 

200 

141 

1969 

98 

110 

108 

46 

1929 

90 

150 

133 

1970 

98 

102 

98 

55 

1930 

84 

108 

108 

1971 

97 

97 

95 

49 

1931 

88 

119 

123 

1972 

98 

93 

90 

46 

1932 

98 

135 

126 

1973 

117 

120 

125 

85 

108 

1933 

93 

119 

119 

1974 

108 

121 

133 

63 

75 

1934 

77 

103 

88 

1975 

133 

114 

132 

49 

78 

1935 

122 

171 

139 

1976 

119 

97 

118 

46 

62 

1936 

91 

115 

92 

1977 

112 

90 

102 

48 

70 

1937 

105 

167 

107 

1978 

107 

86 

114 

43 

64 

1938 

117 

174 

129 

1979 

87 

93 

104 

45 

83 

1939 

95 

140 

99 

I960 

81 

83 

90 

66 

76 

1940 

86 

109 

93 

1981 

72 

66 

78 

52 

75 

1941 

73 

86 

74 

1982 

104 

87 

99 

40 

72 

1942 

90 

127 

106 

1983 

114 

103 

98 

40 

76 

1943 

79 

112 

106 

1984 

100 

107 

112 

69 

106 

1944 

100 

101 

91 

1985 

133 

124 

128 

66 

88 

1945 

116 

134 

123 

1986 

121 

116 

118 

62 

92 

1946 

103 

150 

122 

1987 

107 

95 

94 

57 

80 

1947 

97 

150 

134 

1988 

75 

81 

78 

85 

107 

1948 

107 

129 

112 

1989 

74 

102 

lo  i 

S3 

92 

1949 

95 

141 

118 

1990 

74 

96 

107 

50 

106 

1950 

95 

122 

115 

1991 

83 

103 

120 

60 

112 

1951 

118 

161 

138 

101 


